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ABSTRACT 


This  report  discusses  the  theory,  design,  implementation  and  testing  of  a 
personal  computer-based  Multi-Frequency  Modulation  (MFM)  packet  commu¬ 
nications  system.  Transmitter/receiver  programs  provide  software  drivers  for 
D/A  and  A/D  boards  and  perform  symbol  encoding,  modulating,  demodulating 
and  decoding.  The  design  and  construction  of  a  polarity  coincidence  correlator  for 
receiver  packet  synchronization  is  presented.  Experimental  results  show  that  the 
implemented  MFM  communication  system  conforms  to  theoretical  analysis  with 
acceptable  bit  error.  Results  also  show  that  MFM  can  be  uniquely  adapted  to  a 
specific  channel. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research 
may  not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been 
made,  within  the  time  available,  to  ensure  that  the  programs  are  free  of  compu¬ 
tational  and  logic  errors,  they  cannot  be  considered  validated.  Any  application 
of  these  programs  without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


Reliability  and  flexibility  are  two  fundamental  advantages  of  digital  commu¬ 
nications.  They  are  achieved  through  the  speed  and  processing  power  of  digital 
integrated  circuits.  Multi-Frequency  Modulation  (MFM)  uniquely  harnesses  this 
power  found  in  the  modern  personal  computer  to  encode,  modulate,  demodulate 
and  decode  digitally  formatted  information,  for  sample  rates  currently  up  to 
100,000  samples  per  second.  With  computer-to-computer  communication  links 
implemented  in  a  variety  of  mediums,  such  as  wire,  optical  fiber,  radio  frequency, 
or  acoustic,  MFM  readily  adapts  to  a  given  medium  and  can  emulate  most  ex¬ 
isting  signal  modulation  formats. 

The  focus  of  this  thesis  is  the  implementation  of  a  MFM  packet  communi¬ 
cations  system  using  industry  standard  personal  computers  (PC).  Because  of  the 
limited  processing  speed  of  the  PC,  modulation  and  demodulation  of  the  packet 
is  performed  in  longer  than  real  time.  During  the  course  of  this  project,  trans¬ 
mitter  and  receiver  software  were  developed  to  process  data  for  establishing  a  link 
between  PCs  using  D  A  and  A/D  data  acquisition  boards.  A  correlator  was  de¬ 
signed  and  constructed  to  synchronize  the  receiver  to  the  transmitted  packets  us¬ 
ing  a  synchronization  baud  at  the  beginning  of  the  packets. 

Analysis  of  the  system  phase  response  led  to  signal  modifications  which  re¬ 
sulted  in  improved  signal  quality.  Signal-to-noise  performance  was  evaluated  on 
the  system  for  various  packet  constructions.  The  scope  of  the  report  encompasses 
the  theory  of  MFM,  a  description  of  its  specific  development  and  implementa¬ 
tion,  and  a  summary  of  the  evaluation  of  the  system's  performance. 
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II.  THEORY  OF  MULTI-FREQUENCY  MODULATION 


A.  PACKET  CONSTRUCTION 

The  MFM  signal  set  consists  of  "packets"  of  multiple  tones  which  are  ampli¬ 
tude  and/or  phase  modulated.  These  tones  are  present  simultaneously  during  a 
subinterval  of  the  packet  known  as  a  baud.  Packets  can  be  located  arbitrarily  in 
the  frequency  spectrum  and  time  as  seen  in  Figure  1. 

The  follow  definitions  are  used  in  MFM  [Ref.  1:  pp.  5-6]: 

T :  Packet  length  in  seconds 
A 7':  Baud  length  in  seconds 

Ay  Baud  length  in  number  of  samples 
<f),k:  Symbol  set.  </>,*  is  the  phase  of  the  kth  tone  in  the  Ith  baud 

L :  Number  of  baud  per  packet 
At:  Time  between  samples  in  seconds 

fx  =  1/A t:  Sampling  or  clock  frequency  for  D/A  and  A/D  conversion  in  Hz 
A f  =  1 1  AT:  Frequency  spacing  between  tones 

K:  Number  of  MFM  tones 

To  ensure  the  packet  can  be  uniquely  represented  by  the  kx  samples,  the 
Sampling  Theorem  states  that  fx,  the  sampling  frequency,  must  be  greater  then 
twice  the  highest  frequency  in  the  signal  set  [Ref.  2:  pp.  46-56].  Since  fx  =  kxAf 
and  Af  is  fixed  by  AT,  the  highest  tone  is  kJ2  -1.  Therefore,  the  signal  set  can 
consist  of  an  arbitrary  selection  of  tones  from  dc  to  fJ2  -  Af. 

The  analog  representation  for  the  MFM  signal  during  the  Ith  baud  is 

kxn 

x  ft)  =  Y_,  Aik  cos(2ttAA//  +  (/-  \)AT<t<lAT.  ( 1 ) 

k= 0 

The  first  baud  begins  at  t  =  0  and  the  last  baud,  L,  ends  at  T  =  LAT  .  Sampling 
*,(/)  at  t  =  nAt,  where  n  is  discrete  time,  and  substituting  Af-ljAT  and 
At  =  AT/kx,  produces  a  sampled  output  sequence  of 


I< - L - H 


K -  T  - T - >1 

Figure  1.  MFM  Signal  Packet  (after  Ref.  1:  p.  3.) 


kxl  2 

*/(")  =  Yj  A‘k  C0S 

k= 0 

a  digital  signal  sequence  of  length  kx  samples.  The  kx  point  Discrete  Fourier 
Transform  (DFT)  of  xt(n)  is 

kxl2 

X,{k')  =  Yj\kxAik{^,kW  -Q  +  e~j*H{W  -  {kx  -  k))},  0  <;  k'  <  kx  -1.(3) 
*=o  z 


2nkn 


+  < fiik  )j  0  <n<,k 


x> 


(2) 


From  (3),  it  can  be  seen  that  the  upper  half  of  X,(k')  is  the  complex  conjugate 
image  of  the  lower  half,  which  is  consistent  with  the  symmetry  property  of  a  real 
sequence  [Ref.  2:  p.  402]. 
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B.  SIGNAL  GENERATION  AND  DEMODULATION 

In  the  past,  generation  of  x,(/)  was  attempted  by  constructing  multiple  phase 
lock  loops,  which  were  harmonically  related.  In  December  1985,  LCDR  Deborah 
E.  DeFrank,  then  at  NPS,  designed  a  coherent  multifrequency  synthesizer  as  a 
First  step  in  producing  x,(t)  [Ref.  3].  This  method  proved  to  be  difficult  to  im¬ 
plement  due  to  the  phase  jitter  in  the  phase  lock  loops  and  changes  required  in 
hardware  in  order  to  change  the  harmonic  frequencies  used  in  the  tone  set. 

However,  using  the  properties  of  DFT's,  a  baud  of  x,(n)  can  be  generated  by 
the  host  transmit  computer  by  loading  he  first  half  of  a  complex  valued  array 
(0  to  kJ2  -1)  with  the  magnitude  and  phase  of  the  tones  to  be  included.  To  en¬ 
sure  x t(n)  is  real,  the  upper  half,  (kj 2)  to  ( kx  -  1).  is  loaded  with  the  complex 
conjugate  of  the  values  in  the  first  half  of  the  array  at  the  image  harmonics.  An 
example  of  the  DFT  of  a  real  x(n)  having  a  period  of  ^=16  and  3  tones  is  shown 
in  Figure  2.  The  Inverse  Discrete  Fourier  Transform  (IDFT)  generates  the  real 
discrete  sequence,  x,(n).  which  is  clocked  out  of  the  computer  thru  a  D;A  con¬ 
verter,  at  fx  samples  per  second.  A  signal  packet  is  generated  by  L  repetitions  of 
the  above  process. 

Demodulation  of  MFM  is  the  inverse  of  the  signal  generation  process.  The 
analog  signal,  x(t).  is  sampled  at  fx  samples  per  second  and  converted  to  digital 
format  with  an  A/D  converter.  The  sampled  values  are  loaded  into  the  real 
components  of  a  kx  point  complex  array,  with  the  imaginary  component  set  to 
zero.  A  DFT  is  computed  of  the  array  to  obtain  the  complex  frequency  repre¬ 
sentation  of  the  sampled  input.  Since  ihe  upper  half  of  the  DFT  is  redundant 
information,  only  the  lower  half  is  retained  for  further  processing. 

C.  PROPERTIES  OF  MFM 

Several  important  properties  of  MFM  signals  are  presented  in  this  section. 

( 1 )  Orthogonality  of  signal  during  a  baud.  In  continuous  time 


f(i/2)/!*Ar 

lo 


k  =  i 
k*i 


(4) 


4 


k 

Re(X(k)) 

Im(X(k)) 

0 
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0 
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0  • 
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XR3 
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Figure  2.  Data  Structure  (after  Ref.  1:  p.  9.) 
and  in  discrete  time 


n=  1<X  -1 

xk("Mf  0  = 

n=0 


(1/2  )Alkx 

o 


k  =  i 
k*i 


(5) 


The  derivation  of  the  orthogonality  of  harmonically  related  signals  is  well  known 
[Ref.  4:  pp.  152-154].  Two  advantages  of  an  orthogonal  signal  set  are:  l)lhe 
noise,  assuming  additive  white  Gaussian  noise(AWGN),  affects  each  transmitted 
tone  independently,  thus  simplifying  statistical  computations,  and  2)no  tone  in  a 
baud  interferes  with  any  of  the  other  tones.  This  is  shown  by  representing  each 
tone  in  (1),  neglecting  phase  and  magnitude,  as 

sk{t )  =  cos{2nkAft)  rect(t/AT)  (6) 

where 

I 

rect(f/A7)  = 


0  £  t  <,  AT 
otherwise 


(7) 
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The  Fourier  transform  of  sk(t)  is 

Sk{f)  =  A7sine(/'—  kAf)AT.  (8) 

The  spectrum  of  three  adjacent  tones,  plotted  in  Figure  3,  shows  the  peak  of  each 
tone  coinciding  with  zero  crossings  of  the  spectrum  of  all  other,  thus  minimizing 
their  mutual  interference. 

(2)  Autocorrelation  function  (acf).  Assuming  x(n)  is  periodic, 
the  circular  acf  of  (2)  is 

rxip)  =  x(n)x(n®p),  0<p<kx-\  (9) 

«= o 

where  ©  is  a  left  circular  shift.  For  a  white  bandlimited  sequence  with  an  even 
number  of  harmonics,  the  acf  is  given  by 

rx(p)  =  (1 12)  A  \  cos(2*/cop/*x)  0  °) 

where  k0  =  (k{  +  k2)l 2,  the  midband  harmonic,  and  K  =  /c,  -  k2  +  1,  the  number 
of  tones  in  the  baud.  Figure  4  shows  the  acf  of  a  bandpass  sequence  with 
kx  =  256,  /c,  =  68  and  k2  =  S3(£  =  16,  =  75.5).  These  signal  parameters  con¬ 

form  to  those  of  the  synchronization  baud,  which  will  be  discussed  later.  Note  the 
peak  of  the  acf  occurs  at  p{ 0),  this  feature  is  fundamental  to  synchronizing  the 
receiver  to  the  incoming  packet. 

13)  Matched  filter.  It  is  known  that  a  matched  filter  maximizes 
the  signal-to-noise  ratio  for  additive  white  noise  [Ref.  4:  pp.  88-89].  The  DFT 
of  x{n)  at  the  frequency  k  is  identical  to  the  output  of  a  filter  matched  to  xA(/i). 
This  is  due  to  the  orthogonality  property  of  the  MFM  signal  and  the  linearity  of 
the  DFT  [Ref.  1:  pp.  17-21]. 

D.  MODULATION 

Modulation  is  the  process  of  encoding  the  source  information  onto  a  bandpass  signal 
with  a  canier  frequency  f.  This  bandpass  signal  is  caiied  the  modulated  signal  s(/), 
and  the  baseband  source  signal  is  called  the  modulating  signal  m(t)  [Ref.  5:  p.  204]. 
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68.0  69.0  70.0  71.0  72.0  73.0 

TONES 


Figure  3.  Orthogonal  tone  spacing. 

The  MFM  signal  can  uniquely  accomplish  modulation  in  a  number  of  ways. 
The  signal  x,(/i),  as  defined  in  (2),  can  be  modulated  in  amplitude,  frequency,  and 
phase,  by  translating  the  message  into  changes  in  A,k,  k  and  <£,*  respectively.  Note, 
any  combination  of  the  modulation  types  is  also  possible,  such  as  changing  am¬ 
plitude  and  phase  to  produce  quadrature  amplitude  modulation  (QAM).  The 
signal,  j(r),  is  called  a  bandpass  signal.  However,  MFM  can  be  bandpass  or 
baseband  and  through  multiplication  with  a  carrier  frequency  x,(n)  can  be  trans¬ 
lated  to  any  frequency  band  desired.  The  signal  sets  considered  in  this  thesis  are 
bandpass  and  modulated  using  quadrature  phase  shift  keying  (QPSK)  and  dif¬ 
ferential  quadrature  phase  shift  keying  (DQPSK). 
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Figure  4.  ACF  of  a  white  bandpass  sequence  (after  Ref.  I:  p.  16.) 


( 1 )  QPSK.  Conventional  QPSK  converts  a  digital  input  into 
four  modulation  voltage  levels(symbols)  to  determine  the  phase  of  the  transmitter 
output.  To  minimize  the  probability  of  symbol  error  the  phases  are  spaced  at 
multiples  of  n/2  .  A  plot  of  the  complex  envelope  of  one  tone  of  x,(n)  is  shown  in 
Figure  5.  The  angle  <j)  can  take  on  the  values  of  ±  n/4  and  ±  3tt/4. 

Encoding  MFM  with  QPSK  is  accomplished  by  loading  the 
complex  frequency  domain  array  with  the  appropriate  phase  information.  For 
example,  if  the  digital  input  is  '0110../.  The  first  symbol,  '01',  would  be  loaded 
into  the  frequency  bin,  k,  as  Re[A'(&)]  =  and  lm[A'(/c)]  =  Aj-JT,  where 

Ak  is  the  amplitude  of  tone  k.  Likewise,  the  second  symbol  TO'  would  be  loaded 
into  the  second  bin,  k  +  l,  as  Re[A’(&  + 1)]  =  -Ak+l/y/2  and 
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Imaginary 


Real 

(in  phase) 


Figure  5.  QPSK  signal  constellation. 


Im[A'(/c  +  1)]  =  Ak+llJ 2  .  This  would  continue  until  all  the  tones  in  the  baud  were 
filled  or  there  were  no  more  symbols.  As  mentioned  earlier,  the  discrete  time  do¬ 
main  signal  is  produced  by  taking  the  IDFT  of  the  complex  frequency  array. 
Each  successive  baud  is  encoded  similarly. 

Decoding  QPSK  directly  into  bits  is  accomplished  easily  as  fol¬ 
lows.  Assuming  a  coherent  receiver,  decoding  requires  evaluating  the  polarity  of 
the  real  and  imaginary  components  of  each  frequency  bin.  Notice  in  Figure  5 
that  the  symbol  mapping  uses  Gray  encoding.  This  reduces  the  probability  of 
bit  error  because  errors  caused  by  AWGN  are  likely  to  occur  when  the  adjacent 
symbol  is  selected  for  the  transmitted  symbol;  thus,  the  symbol  error  will  contain 
only  one  bit  error.  Gray  encoding  also  allows  decoding  straight  into  bits,  with  the 
right  bit  determined  by  the  polarity  of  the  real  component  and  the  left  bit  by  the 
polarity  of  the  imaginary  component.  The  digital  signal  is  obtained  by  succes¬ 
sively  decoding  each  tone  of  each  baud. 

(2)  DQPSK.  QPSK  requires  strict  phase  coherent  regeneration 
of  the  sampling  frequency  to  eliminate  phase  ambiguity.  This  results  in  a  complex 
synchronization  design  dr  distribution  of  the  sampling  clock  frequency  to  each 
receiver.  DQPSK  provides  a  practical  solution  to  the  phase  uncertainty  problem 
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but  at  a  cost  of  approximately  2.3  dB  more  in  SNR  in  order  to  obtain  the  same 
probability  of  bit  error  as  coherent  QPSK  [Ref.  6:  p.  442],  Encoding  MFM  with 
DQPSK  is  similar  to  QPSK;  however,  DQPSK  translates  the  original  symbol  set 
into  a  second  "differential"  symbol  set,  which  is  then  encoded  as  QPSK.  Trans¬ 
lated  symbols  are  generated  based  on  the  input  symbol  and  the  previous  trans¬ 
lated  symbol.  Figure  6  shows  this  translation.  Notice  that,  regardless  of  the 
previous  symbol,  an  input  of  W  generates  a  new  symbol  in  the  same  quadrant 
as  the  previous  symbol.  ’An  input  of  '01'  rotates  the  new  symbol  +  n/2  radians 
from  the  previous,  Tl'  rotates  n  radians,  and  TO'  rotates  -  n/2.  In  the  receiver, 
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decoding  is  performed  by  determining  the  phase  difference  between  successive 
pairs  of  tones,  using  complex  arithmetic.  The  phase  difference,  A 4>,  is  found  from 

^  =  =  e^'  (11) 

and  in  the  absence  of  noise  will  be  0,  ±  rc/2,  or  n  radians.  To  realign  the  complex 
signal  to  the  original  constellations,  e^6  is  rotated  by  +  n/4  radians  and  the  in- 
phase  and  quadrature  components  are  decoded  as  with  QPSK. 
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III.  SYSTEM  DEVELOPMENT 


Multi-Frequency  Modulation  joins  the  basic  methods  of  digital  communi¬ 
cations  theory  and  signal  processing.  The  theory  behind  MFM  is  not  new  [Ref. 
7],  however,  the  ability  to  implement  it  inexpensively  is  linked  to  the  recent  de¬ 
velopment  of  inexpensive  DSP  chips  for  personal  computers.  The  advent  of 
packet  switching  in  data  communications  also  makes  MFM  a  preferred  choice  for 
the  MODEM  because  of  its  "packet"  like  format. 

A.  SYSTEM  DESCRIPTION 

An  MFM  communications  system  is  shown  in  Figure  7.  Information,  de¬ 
noted  by  m{n ),  is  modulated  into  a  frequency  band  that  will  propagate  over  the 
available  channel.  The  receiver  converts  the  noise-corrupted  signal,  ;*(r),  to  an 
estimate  of  the  source  information,  m(n).  Theoretically,  m(n)  is  any  signal  that 
can  be  represented  in  a  digital  format.  To  be  practically  useful,  the  conversion 
process  must  be  real  time.  However,  due  to  signal  conversion  speeds  available  for 
this  research  project,  the  input  information  had  to  be  restricted  to  stored  data 
files.  The  ideal  transmitter  output,  x(/),  is  described  for  the  l,h  baud  by  (1). 


n(t) 


Figure  7.  MFM  Communication  System. 


B.  BLOCK  DESCRIPTION 
1.  Transmitter 

The  transmitter  is  subdivided  into  functional  blocks,  shown  in  Figure  8. 

I 

As  mentioned  above,  m(n)  has  been  digitized  and  stored  in  a  file.  Therefore,  the 
input  is  a  string  of  binary  digits.  Processing  is  performed  on  a  baud-by-baud 
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Figure  8.  Transmitter  Functional  Block  Diagram. 


basis,  until  the  end  of  the  data  file,  or  maximum  packet  length  has  been  obtained, 
whichever  comes  first. 

The  encoder  converts  input  symbols  into  complex  values  stored  in  the 
frequency  domain  array.  The  value  of  the  symbols  depends  on  the  type  of  mod¬ 
ulation.  For  exampie,  QPSK  symbols  are  two  bits  long  and  are  encoded  as  pre¬ 
viously  discussed.  The  discrete  signal,  produced  by  computing  the  IDFT  of  the 
complex-valued  frequency  domain  array,  is  loaded  into  the  packet  storage  area. 
The  controller  determines  the  parameters  of  the  packet  based  on  the  modulation 
type,  baud,  and  message  size;  it  then  sequences  the  input  data  through  the 
transmitter  one  baud  at  a  time.  Once  the  message  has  been  processed,  the  entire 
stored  digital  signal,  jc(n),  is  transferred  out  at  the  selected  rate,  fx  samples  per 
second,  through  a  D/'A  converter.  Depending  on  the  channel's  frequency  re¬ 
sponse,  filtering  of  the  output  may  be  desired  to  remove  the  high  frequency  am¬ 
plitude  discontinuities  introduced  by  the  D/A  converter. 

2.  Receiver 

The  receiver,  shown  in  Figure  9,  demodulates  the  MFM  signal  by  re¬ 
versing  the  transmitter  process.  Data  acquisition  is  the  process  that  samples  the 
analog  signal  at  fx  samples  per  second  and  converts  it  to  digital  format  with  an 
A/D  converter.  Though  not  shown  in  Figure  9,  filtering  of  the  input  is  recom¬ 
mended  to  bandlimit  input  noise  and  to  reject  out-of-band  interference.  The 
converted  data  sequence  is  stored  until  all  data  is  acquired.  Again,  this  is  due  to 
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time  limitations  in  the  signal  processing  algorithms  available  in  the  system  used 
in  this  research  task.  The  stored  real  values  are  accessed  one  baud  at  a  time  to 
perform  a  k,  point  DFT.  The  first  half  of  the  resulting  complex  values  are  de¬ 
coded  to  obtain  the  amplitude  and  phase  modulation  information.  As  in  the 
transmitter,  the  controller  sequences  the  data  through  the  system. 


Figure  9.  Receiver  Functional  Block  Diagram. 

Notice  the  delay  from  input  to  output  is  the  time  it  takes  to  acquire  the 
entire  packet  and  process  the  first  baud.  In  contrast,  real  time  processing  must 
complete  the  kx  point  DFT  in  the  same  or  less  time  than  it  takes  to  fill  up  one 
baud  length  buffer.  Processing  can  alternate  between  buffers,  and  data  will  flow 
continuously  through  the  system  at  the  sampling  rate.  The  storage  requirements 
and  processing  time  usings  =  61440  Hz,  L  =30,  and  kx  =  1024  samples  are  given 
in  Table  1  for  a  system  with  a  real  time  and  non-real  time  DFT  capability  and 
for  the  non-real  time  system  used  in  this  research.  The  AT  implementation 
processing  time  highlights  the  need  for  a  high  speed  vector  processor  to  perform 
the  DFT  computations. 

3.  Synchronization 

In  digital  communications,  various  degrees  of  synchronization  are  re¬ 
quired.  Typical  degrees  are: 

♦  Carrier  synchronization, 

*  Bit  or  symbol  synchronization, 

•  Word  or  frame  synchronization. 
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Tabic  1.  STORAGE  REQUIREMENTS  AND  PROCESSING  TIME 


Real  time 

Non  real  time 

AT  implementa¬ 
tion 

Storage  required 
(#  of  samples) 

2048 

30720 

30720 

Time  required  for 
processing  1024 
samples  (msec) 

<:  16.7 

>  16.7 

10  x  103 

These  have  slightly  different  meanings  and  forms  depending  on  the  system.  For 
the  MFM  system,  carrier  and  frame  synchronization  may  be  required. 

Carrier  synchronization  is  required  if  the  modulation  scheme  is  coherent, 
such  as  QPSK.  This  means  generating  an  fx  at  the  receiver  in  frequency  and  phase 
coherence  with  the  transmitted  fx.  Though  the  packet  does  not  contain  the  fx 
harmonic,  every  tone  is  harmonically  related  to  it.  The  packet  can  be  constructed 
to  have  a  pilot  tone  separated  from  the  modulation  tones  from  which  fx  can  be 
derived.  For  example,  the  modulation  tones  could  range  from  k,  to  kH  ,  giving  a 
band  of  j  +  1  consecutive  tones.  The  pilot  tone,  kp,  would  then  be  placed  several 
tones  away  to  ensure  minimum  band  interference  from  the  message  during  its 
extraction. 

As  mentioned  previously,  MFM,  as  implemented,  is  in  the  form  of  a 
transmission  packet.  To  acquire  the  packet  the  receiver  must  know  when  to  start 
sampling.  This  is  accomplished  by  frame  synchronization.  Typically,  unique 
words  are  inserted  to  mark  the  start  of  each  frame.  In  MFM,  the  unique  word  is 
called  a  "synchronization  baud",  and  it  is  added  at  the  beginning  of  each  packet. 
This  baud  is  generated  similarly  to  other  bauds,  except  the  tones  and  phases  are 
predetermined.  Acquisition  of  the  received  signal  starts  after  successful  detection 
of  the  synchronization  baud.  [Ref.  5:  pp.  511-512,  8:  pp.  293-295] 
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IV.  SYSTEM  IMPLEMENTATION 


In  the  previous  chapters  the  tools  and  framework  for  developing  MFM 
packet  signals  were  established.  One  actual  realization  of  this  form  of  communi¬ 
cations  will  now  be  presented.  Though  the  IBM  Personal  Computer  (PC)  was 
used,  MFM  can  be  implemented  in  a  variety  of  ways.  The  PC,  however,  is  ide¬ 
ally  suited  for  MFM  implementation  for  the  following  reasons: 

•  Increased  I/O  channel  maximum  throughput  rate  of  approximately 
100  KHz  using  direct  memory  access  (DMA). 

•  Digital  signal  processing  allows  encoding,  decoding,  modulation,  demodu¬ 
lation,  and  channel  equalization. 

•  Signal  processing  algorithms  are  available  in  high  level  languages. 

•  Packet  construction  is  easily  modified  to  conform  to  various  channels. 

•  External  hardware  is  easily  interfaced. 

•  Cost  is  low. 

Successful  data  acquisition  and  decoding  of  (1)  was  the  first  and  most  im¬ 
portant  goal  of  this  thesis.  Other  goals  affecting  software  and  hardware  develop¬ 
ment  were  the  following: 

•  Maximize  packet  size. 

•  Develop  synchronization  circuitry. 

•  Transmit, receive  ASCII  files. 

•  Develop  software  for  testability  and  flexibility. 

•  Develop  software  for  statistical  testing. 

A.  SIGNAL  PARAMETERS 

Software  and  hardware  implementations  can  be  easily  modified  to  receive  any 
packet  construction.  However,  all  software  conforms  to  the  signal  parameters  in 
Appendix  A,  established  between  NPS  and  NOSC  for  an  acoustic  application. 
This  is  a  bandpass  signal  in  the  band  from  16-20  KHz  using  a  clock  rate  of 
61440  Hz.  All  baud  lengths  are  powers  of  two,  allowing  utilization  of 
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hardware/software  Fast  Fourier  Transform  (FFT)  algorithms  to  speed  frequency 
transformation  processing. 

Various  baud  sizes  are  available  to  give  greater  flexibility  in  adapting  the 
packet  to  test  specific  channel  parameters.  In  DQPSK,  where  the  information  is 
represented  by  the  phase  difference  between  adjacent  tones,  channel  phase  dis¬ 
tortion  affects  the  shorter  baud  more,  due  to  their  larger  A/.  The  maximum 
tone-to-tone  phase  error  introduced  by  a  linear  phase  channel  is  shown  in 
Table  2. 


Table  2.  PHASE  ERROR  FROM  A  LINEAR  CHANNEL 


256 

512 

1024 

2048 

4096 

A/ (Hz) 

240 

120 

60 

30 

15 

Phase  error 
(degrees) 

25.6 

12.8 

3.2 

1.6 

Obviously  the  longer  baud  would  be  preferred  when  differential  coding  is  between 
adjacent  tones.  However,  if  the  channel  introduces  time-related  distortion,  like 
noise  bursts,  electrical  glitches,  or  propagation  fluctuations,  a  short  duration  baud 
is  desired,  and  differential  coding  should  be  between  the  same  tones  on  adjacent 
bauds. 

B.  SOFTWARE 
1.  Transmitter 

Generation  of  the  MFM  signal  has  been  accomplished  in  previous  re¬ 
search  at  NPS.  Utilizing  previous  hardware  and  core  software,  this  project  ex¬ 
panded  the  transmitter  software  to  provide  encoding  of  16-QAM  and  QPSK, 
encoded  data  files  using  DQPSK,  constructed  maximum  size  packets,  generated 
a  synchronization  baud,  and  provided  greater  flexibility  in  initializing  the  DMA. 
Actual  transmitter  hardware  will  not  be  discussed;  it  has  been  the  subject  of 
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previous  research  [Ref.  9].  Two  new  transmitter  programs,  TRANSMIT  and 
XMITMES,  will  be  discussed. 
a.  TRANSMIT 

The  program  TRANSMIT  provides  maximum  flexibility  in  con¬ 
structing  and  encoding  signal  packets.  It  is  an  excellent  tool  for  basic  system 
testing  and  can  be  used  as  a  training  aid  to  demonstrate  various  digital  modu¬ 
lation  schemes.  A  procedural  flow  diagram  of  TRANSMIT  is  shown  in 
Figure  10.  SeiectBaud  is  a  user  interactive  procedure  for  establishing  the  pa¬ 
rameters  of  a  16-20  KHz  packet  and  for  choosing  the  encoding  scheme  to  be 
utilized.  Parameters  kx  and  L  are  selected.  From  kx,  the  lower  and  upper 
bandlimits,  and  k2,  are  set.  This  generates  a  bandlimited  packet  with 
k2  —  k\  +  1  tones.  The  data  storage  requirement  is  kxL,  representing  the  total 
number  of  samples  in  the  packet  which  are  clocked  out.  The  encoding  scheme 
determines  the  path  for  encoding  the  cpmplex  frequency  domain  array. 


EncodeQPSK  begins  by  displaying  a  four  symbol  QPSK  constellation 
used  as  a  reference  in  selecting  the  symbols  over  the  band.  The  symbols  for  the 
K  tones  in  the  baud  can  be  selected  in  a  variety  ways: 

•  Symbols  for  all  tones  are  randomly  selected  from  a  random  generator. 

•  The  symbol  for  each  lone  is  selected  by  the  user. 

•  Individual  tones  may  be  removed  from  the  band. 


18 


This  last  feature  allows  construction  of  a  baud  with  an  arbitrary  number  of  tones 
within  the  band  determined  by  A,  and  A2.  Selected  symbols  for  the  band  are 
loaded  into  the  complex  frequency  domain  array  with  their  complex  conjugate 
image  frequency.  EncodeQAM  is  functionally  the  same  as  EncodeQPSK,  except 
the  symbols  conform  to  a  16-QAM  constellation  [Ref.  1:  p.27j. 

To  obtain  the  real  discrete  time-domain  sequence  of  the  encoded 
baud,  ComplexFFT  computes  the  inverse  FFT.  ComplexFFT  consumes  the  ma¬ 
jority  of  the  processing  time  in  the  program  due  to  the  complex  arithmetic  oper¬ 
ations  required,  thus  restricting  the  overall  throughput  of  the  system.  For  real 
time  processing,  the  FFT  algorithm  must  be  accomplished  by  a  hardware  signal 
processor. 

Each  value  in  the  time  domain  sequence  is  represented  as  a  real  data 
type,  occupying  six  bytes  of  memory.  ScaleData  converts  these  values  down  to  a 
one  byte  format  acceptable  to  the  D/A  converter  and  places  them  into  a  packet 
storage  buffer.  EncodeQPSK,  ComplexFFT,  and  ScaieData  are  executed  for 
each  baud,  until  all  L  baud  have  been  processed.  To  transmit  the  packet  out  of 
the  computer,  DMAINIT  transfers  samples  at  fx  samples  per  second  over  the 
DMA  channel  to  D;A  converter  [Ref.  9,10]. 
b.  XMITMES 

To  demonstrate  the  suitability  of  MFM  for  transferring  information 
from  a  source  to  a  sink,  the  program  XMITMES  was  written  to  transmit  an 
ASCII  file  encoded  using  DQPSK.  Affixed  to  the  beginning  of  the  packet  is  a 
synchronization  baud.  As  shown  in  Figure  11,  XMITMES  has  an  even  simpler 
structure  then  TRANSMIT,  because  all  tones  in  the  band  are  encoded  using 
DQPSK. 

The  synchronization  baud  is  a  predetermined  sequence  generated  by 
SvncBaud.  This  baud  is  constructed  as  are  all  other  baud,  except  kx  is  fixed  at 
256,  and  tones  68  to  83  are  encoded  with  the  same  random  symbol  pattern  re¬ 
gardless  of  the  packet  construction  or  input  message.  This  synchronization  se¬ 
quence  occupies  the  first  256  values  in  the  packet  buffer  and  therefore  is  the  first 
to  be  clocked  out  of  the  computer. 
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As  messages  can  be  of  various  sizes,  TailorPacket  sets  the  maximum 
number  of  baud  required  for  encoding.  This  is  determined  by  dividing  the  num¬ 
ber  of  characters  in  the  message  by  the  number  of  characters  that  can  be  en¬ 
coded. 

DiffEncode  encodes  the  message  file  into  the  complex  frequency  do¬ 
main  array.  It  reads  one  character  at  a  time;  then  breaks  the  eight  bit  character 
into  four  2-bit  symbols.  The  symbols  are  DQPSK  encoded  and  stored  in  the 
frequency  array.  Once  encoded,  processing  and  signal  output  by  ComplcxFFT, 
ScaleData,  and  DMAIN1T  are  the  same  as  in  TRANSMIT. 

2.  Receiver 

Like  the  transmitter,  the  receiver  requires  hardware  to  interface  the  com¬ 
puter  with  the  channel  and  processing  software  to  demodulate  the  MFM  signal 
packet.  The  channel  interface  is  a  high-speed  data  acquisition  board,  model 
DASH-16F,  manufactured  by  MetraByte  Corporation  [Ref.  11].  Also  at  the  re¬ 
ceiver,  synchronization  circuitry  is  required  to  detect  the  beginning  of  the  packet. 
The  synchronization  circuitry  is  discussed  in  Section  C.  Implementation  of  the 
processing  software  assumes  synchronization  and  a  receiver  fx  the  same  as  at  the 
transmitter.  Two  receiver  programs,  RECEIVE  and  RECMES  will  be  discussed. 
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a.  RECEIVE 


The  program  RECEIVE,  shown  in  Figure  12,  processes  and  displays 
only  the  first  baud  of  a  packet.  It  is  an  excellent  tool  to  experiment  with  the 
DASH-l<3F's  software  drivers  and  it  provides  a  quick  indication  of  the  overall 
system  performance.  Two  procedures  are  available  for  data  input,  GetData  and 
AcquireData.  GetData  reads  in  the  time  domain  sequence  File  generated  in  the 
transmitter,  eliminating  all  data  acquisition  hardware  in  the  transmitter  and  re¬ 
ceiver;  thus  ensuring  perfect  synchronization  and  accurate  evaluation  of  all  en¬ 
coding,  decoding,  and  processing  software.  The  computational  SNR  of  the  system 
is  determined  under  this  condition,  as  the  only  noise  present  is  from  round  off 
and  truncation  errors  introduced  in  processing. 


Figure  12.  RECEIVE  algorithm. 


Analog  data  acquisition  is  performed  by  the  procedure  Acquiredata. 
It  initializes  and  controls  the  DASH-1 6F  using  procedures  written  by  Quinn- 
Curlis  [Ref.  12].  AcquireData  allocates  memory  to  store  the  sampled  values 
transferred  from  the  board  using  the  DMA  controller.  Direct  Memory  Access  is 
the  only  data  transfer  mode  capable  of  transferring  data  to  memory  at  the  re¬ 
quired  fx,  without  disruption  by  other  interrupt  processes  in  the  computer.  Other 
important  initialization  parameters  are  triggering  source  and  the  number  of  sam¬ 
ples  to  be  collected.The  A/D  may  be  triggered  from  two  sources,  a  programmable 
interval  timer  or  an  external  trigger  source.  The  programmable  interval  timer 
divides  either  a  1  MHz  or  10  MHz  clock  to  derive  the  sampling  rate  of  the  trig- 
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ger.  Since  this  method  cannot  produce  an  arbitrary  fx,  external  triggering  is  used. 
After  initialization,  conversions  take  place  on  the  positive  transition  of  every 
trigger  until  the  specified  number  of  samples  have  been  acquired  and  transferred 
to  memory.  Collected  data  is  then  converted  into  a  format  acceptable  for  further 
processing. 

DiffDecode  determines  the  encoded  symbols  by  differentially  decod¬ 
ing  the  complex  frequency  array  transformed  from  the  sampled  data  by 
ComplexFFT.  The  output  of  the  program  is  the  frequency  spectrum  plot  in  Fig¬ 
ure  13.  The  lower  plot  is  the  spectral  response  of  the  kJ2  tones  in  the  baud;  the 
upper  plot  represents  only  tones  from  kx  to  k2 .  The  color  on  the  original  display 
indicates  the  phase  quadrant  of  a  given  tone.  Graphing  the  data  in  this  fashion 
provides  quick  qualitative  analysis  of  the  frequencies  and  their  phase  informa¬ 
tion. 


Figure  13.  Baud  magnitude  spectrum, 

b.  RECMES 

RECMES,  shown  in  Figure  14,  demodulates  the  ASCII  encoded 
transmission  produced  by  XMITMES.  It  differs  from  RECEIVE  in  that  it  can 
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Figure  14.  RECMES  algorithm. 

process  a  multiple  baud  packet.  The  user  interactive  procedure,  PacketSetUp, 
tailors  processing  to  the  expected  receive  packet,  using  the  inputs  kx  and  L . 
AcquireData  samples  and  stores  a  memory  segment  of  data  regardless  of  the 
packet  size,  Lkx.  However,  once  stored,  ConvertDala,  ComplexFFT,  and 
DiffDccode  process  only  L  baud  of  the  data.  DiffDecode  combines  four 
differentially  decoded  symbols  into  one  byte,  representing  the  ordinal  number  of 
an  ASCII  character.  To  reconstruct  the  message,  the  characters  are  transferred 
to  iext  file  MESSAGE.DAT  until  processing  is  complete.  For  convenience, 
ShowMessage  displays  the  recovered  message. 

C.  HARDWARE 

Synchronization  of  DQPSK  MFM  is  obtained  from  a  hardware  correlator 
that  is  external  to  the  host  receiver  computer.  The  128  point  correlator,  illus¬ 
trated  in  Figure  15,  provides  the  data  acquisition  board  with  sampling  triggers 
synchronized  with  respect  to  time  of  arrival  of  the  packet.  Using  only  the  polarity 
of  the  analog  input,  it  functions  as  a  matched  filter  to  the  last  half  of  the  256 
point  synchronization  baud.  This  type  of  correlator  is  referred  to  as  a  polarity 
coincidence  corrclator(PCC). 

The  hard  limiter  used  to  obtain  the  polarity  information  in  the  analog  input 
is  a  fast,  high  precision,  high  gain,  operational  amplifier.  During  positive  and 
negative  portions  of  the  input  the  output  is  +5  Vdc  and  0  Vdc  respectively. 
This  unipolar  signal  is  synchronized  to  the  receiver's  fx  as  it  is  clocked  through  a 

P 

128  point  serial  shift  register.  To  minimize  bit  instability,  due  to  the  shift  register 
input  being  asynchronous  to  fx ,  the  hard  limiter  slew  rate  should  be  as  large  as 
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Figure  15.  Polarity  coincidence  correlator.  > 


possible.  For  example,  during  testing,' .15  bits  out  of  128  were  unstable  using  an 
MCI 747  general  purpose  operational  amplifier  with  a  slew  rate  of  0.7  V/jusec. 
This  was  reduced  to  2  bits  using  the  faster  LM318  with  a  slew  rate- of 
70  V/jusec. 

The  l's  and  0's  buses  store  the  time-reversed  polarity  sequence  of  the  syn¬ 
chronization  baud.  The  voltage  on  the  buses  represent  the  correlation  between  the 
stored  sequence  and  the  sequence  in  the  shift  register.  As  the  sequences  come  into 
alignment,  the  l's  bus  voltage  increases  while  the  0's  decreases.  This  inverse 
property  is  combined  by  a  differential  amplifier  to  give  the  total  correlation. 
When  both  sequences  match,  the  l's  bus  voltage  is  +5  Vdc  and  the  0's  is  0  Vdc, 
giving  a  maximum  differential  voltage  of  5  Vdc  into  the  threshold  detector.  The 
threshold  detector  generates  a  synchronization  trigger  on  detection  of  the  corre¬ 
lation  peak.  A  latch  is  set,  enabling  the  D  flip-flop  to  pass  sampling  triggers  for 
data  acquisition  at  the  packet's  beginning  as  required  for  demodulation. 

During  design  and  testing  the  correlator  progressed  from  a  16  to  128  point 
shift  register.  With  each  shift  register  output  bus  connected  through  a  resistor 
voltage  divider  network,  progressive  testing  was  necessary  to  ensure  power  sup¬ 
ply  current  ratings  were  not  exceeded.  Maximal-length  sequences  were  used  as 
test  correlation  sequences  due  to  their  unique  two-valued  autocorrelation  func¬ 
tions  that  are  easily  determined  by  finite-field  arithmetic  [Ref.  8:  pp.  368-375]. 
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V.  PERFORMANCE  EVALUATION 


Initial  testing  of  the  MFM  system  was  concerned  with  quickly  determining 
the  quality  of  the  demodulated  signals.  Figure  13  provided  an  excellent  repre¬ 
sentation  of  the  demodulated  and  decoded  signal  and  was  used  throughout  sys¬ 
tem  development.  This  plot  uniquely  displays  the  frequency  domain  magnitude 
and  phase.  Since  QPSK  symbols  are  at  least  90  degrees  apart  each  color  (on  the 
original  displays)  represents  a  specific  quadrant  in  the  complex  frequency  plot.  For 
example,  a  green  line  indicates  a  tone  whose  phase  is  in  the  second  quadrant  and 
whose  magnitude  is  represented  by  the  length  of  the  line.  Although  each  color 
represents  a  broad  decision  region  in  phase,  it  is  still  clear  a  phase  shift  of  at  least 
ISO  degrees  occurs  across  the  4  KHz  band.  This  shift  is  primarily  due  to  the  linear 
phase  Butterworth  filter  used  to  simulate  the  system  channel,  and  it  contributed 
to  the  decision  to  develop  a  DQPSK  encoded  signal. 

A.  SYSTEM  FREQUENCY  RESPONSE 

After  development  of  all  hardware  and  software,  testing  revealed  significant 
phase  fluctuations  across  the  frequency  band.  The  system  phase  response,  plotted 
in  Figure  16,  was  determined  by  subtracting  the  transmitted  phase  from  the  re¬ 
ceived  phase  for  each  tone.  Ideally  the  response  should  be  smooth.  With 
DQPSK  encoded  from  tone-to-tone,  phase  differences  introduced  by  the  system 
will  result  in  reduced  tolerance  to  additive  noise.  For  example,  notice  the  phase 
difference  of  0.4  radians  between  tones  296  and  297.  When  noise  is  added,  de¬ 
coding  the  symbol  in  error  will  occur  at  a  lower  noise  level  than  it  would  if  the 
phase  difference  was  zero;  thus  a  higher  SNR  is  required  for  low  error  rate  de¬ 
coding  when  the  channel  introduces  tone-to-tone  phase  fluctuations.  The  tone- 
to-tone  phase  fluctuations  were  substantially  reduced,  as  shown  in  Figure  17,  by 
decreasing  the  magnitude  of  the  encoded  tones  in  the  transmitter  which  moved 
the  operating  range  of  the  D/A  converter  to  a  more  linear  region.  Further  anal¬ 
ysis  of  the  effect  of  the  D/A  and  A/D  converters  on  phase  shift  is  recommended 
in  order  to  be  able  to  specify  them  properly  in  future  designs. 
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Figure  16.  Initial  system  response. 

A  second  important  feature  of  the  system  phase  response  is  its  average  slope 
across  the  band.  The  phase  slope  is  affected  by  synchronizatio  timing  as  shown 
in  Figure  18.  When  the  sequence  through  the  correlator  was  delayed  2  bits,  or 
2/61440  =  32.55  usees,  relative  to  the  analog  input  to  the  receiver,  thus  delaying 
data  acquisition  accordingly,  the  phase  slope  changes  as  predicted  by  the  time 
delay  Fourier  transform  theorem; 

x(t-Td)  XWe-^.  (12) 

This  shows  the  phase  slope  in  radians  for  a  signal  delayed  by  Td  seconds,  is  27t7’rf 
radians/'Hz.  For  ’example,  at  /=  Afk,  =  (60)(300)  =  18000  Hz  and 
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Figure  17.  Improved  system  response. 

Td=  32.55  fi secs,  the  phase  shift  is  3.7  radians.  This  compares  closely  to  the  ac¬ 
tual  phase  shift  of  tone  3G0  in  Figure  17  and  Figure  18. 

In  summary,  system  frequency  response  analysis  identified  phase  shifts  in¬ 
troduced  by  the  system  hardware  and  by  the  synchronization  timing.  The  mag¬ 
nitude  of  the  encoded  tones  and  the  synchronization  sequence  delay  were  selected 
experimentally  in  order  to  minimize  tone-to-tone  phase  fluctuations  and  provide 
the  comparatively  flat  response  shown  in  Figure  18  for  all  subsequent  testing. 
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B.  SNR  PERFORMANCE 

The  performance  of  the  system  when  corrupted  by  additive  noise  is  deter¬ 
mined  by  its  output  signal-to-noise  ratio  (SNR).  For  QPSK-encoded  MFM  the 
mean  of  the  real  and  imaginary  parts  of  the  2 K  coefficients  of  the  DFT  represent 
the  received  signal  amplitude,  and  their  variance  represents  the  noise  power.  The 
output  SNR  is  defined  as  the  ratio  of  the  square  of  the  mean  to  the  variance  of 
each  for  these  2K  coefficients  [Ref.  1:  p.  25], 

Figure  19  shows  the  system  SNR,  which  is  the  output  SNR  when  there  is  no 
additive  noise  in  the  channel,  versus  the  frequency  spacing,  A/,  of  the  five  baud 
types.  As  would  be  expected,  DQPSK  performs  better  with  a  smaller  A/because 
the  phase  difference  between  adjacent  tones  is  smaller  when  A/is  smaller. 
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In  Figure  20,  the  SNR  out  of  the  receiver  is  shown  versus  input  SNR  which 
was  varied  by  adding  noise  to  the  analog  signal.  Theoretically,  the  output  SNR 
should  equal  the  narrowband  input  SNR  in  QPSK  [Ref.  1:  pp.  24-25].  The  re¬ 
duction  in  output  SNR  at  higher  input  SNRs  is  due  to  each  baud  approaching 
its  maximum  output  SNR  as  set  by  system  noise  and  shown  in  Figure  19. 
Table  3  lists  the  bit  errors  at  various  SNR  input  levels.  There  were  approxi¬ 
mately  2500  bits  transmitted  for  each  entry.  This  clearly  indicates  the  worsening 
effect  of  the  channel  on  bauds  with  greater  tone  spacing  A f  at  high  SNRs. 

A  direct  comparison  of  the  MFM  signal  with  commercial  modems  is  difficult 
due  to  the  wide  variety  of  techniques  that  are  used  to  achieve  high-speed  data 
transfer.  However,  with  a  bandwidth  efficiency  of  2  bits/s/Hz,  a  DQPSK  MFM 
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Figure  20.  SNR  performance. 


signal  can  achieve  a  throughput  of  6000  bits  per  second  on  a  standard  switched 
telephone  line  having  a  bandwidth  of  3  KHz.  Though  not  implemented  in  this 
project,  several  techniques  could  be  utilized  to  increase  the  bandwidth  efficiency, 
such  as  increasing  the  constellation  size  and/or  using  data  compression  algorithms 
as  in  commercial  modems  [Ref.  13,  14].  Using  these  techniques  MFM  could 
achieve  a  bandwidth  efficiency  of  4-8  Bits/s/Hz,  which  compares  favorably  to 
current  high  speed  modems. 
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Table  3.  BIT  ERRORS  IN  2500  BITS  TRANSMITTED  VS  BAUD  TYPE  AND 
SNR. 


A f 

kx 

SNR  (dB) 

5 

10 

15 

20 

240 

256 

165 

35 

18 

10 

120 

512 

154 

16 

3 

0 

60 

1024 

148 

17 

0 

0 

30 

2048 

140 

9 

0 

0 

15 

4096 

127 

10 

0 

0 
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VI.  CONCLUSIONS  AND  FURTHER  STUDY 


MFM  is  well  suited  to  the  signal  processing  environment  of  the  personal 
computer.  As  implemented  in  this  thesis,  the  packet  signal  provides  a  high  speed, 
low  bit  error  data  transfer  link  between  two  industry  standard  computers.  Syn¬ 
chronization  hardware  allows  the  link  to  use  asynchronous  data. 

The  theory  and  properties  of  MFM  have  been  discussed  for  a  white  bandpass 
packet  signal.  Packet  construction  and  various  modulation  formats  have  been 
shown  to  be  easily  adapted  to  a  given  channel.  The  system  was  developed  and 
implemented  for  testability  and  to  be  easily  modified  to  accommodate  a  variety 
of  applications.  However,  the  data  storage  requirements  and  signal  conversion 
speeds  were  unique  to  this  project  due  to  available  hardware.  As  implemented, 
this  system  would  require  a  multiple  packet  transmission  to  achieve  practical  op¬ 
eration. 

Analysis  of  the  system's  phase  response  led  to  adjustments  in  the  encoded 
tone  magnitudes  and  synchronization  timing.  This  substantially  improved  the 
system's  phase  response.  SNR  test  results  indicate  a  smaller  A/baud  has  superior 
performance  in  a  linear  phase  channel.  For  A/ less  than  60  Hz,  bit  errors  are 
acceptably  low  for  input  SNRs  greater  than  approximately  15  dB.  Exhaustive 
testing  is  required  to  compare  bit  error  performance  with  conditional  digital 
modulation  schemes. 

Areas  of  further  study  should  focus  toward  increasing  the  system  bandwidth 
efficiency  and  implementation  of  hardware  signal  processing  for  real-time  opera¬ 
tion.  In  fact,  current  research  at  NPS  is  directed  at  increasing  the  throughput 
speed  of  the  system  described.  Improved  bandwidth  efficiency  can  be  obtained 
through  data  compression  algorithms  and  developing  receiver  software  to.  decode 
the  16-QAM  encoded  signal  generated  by  TRANSMIT.  Finally  to  provide  a 
greater  flexibility  in  adapting  MFM  to  channels  with  significant  uncompensated 
phase  distortion,  the  option  should  be  developed  to  encode  DQPSK  on  a  baud- 
to-baud  basis. 
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APPENDIX  A.  DESIGN  PARAMETERS 


Table  4.  DESIGN  PARAMETERS  FOR  A  1/15TH  SECOND  SIGNAL  PACKET 
IN  A  16-20KHZ  BANDPASS  CHANNEL. 


Baud  length  (sec) 

AT 

1/240 

1/120 

1/60 

1/30 

1/15 

No.  of  Bauds 

L 

16 

8 

4 

2 

1 

Tone  spacing 

A/ 

240 

120 

60 

30 

15 

Lowest  Harmonic 

68 

135 

269 

537 

1073 

Lowest  tone  freq 

/. 

16320 

16200 

16140 

16110 

16095 

Highest  Harmonic 

83 

166 

332 

664 

1328 

Highest  tone  freq 

h 

19920 

19920 

19920 

19920 

19920 

Samples  per  Baud 

kx 

256 

512 

1024 

2048 

4096 

Sampling  freq 

fx 

61440 

61440 

61440 

61440 

61440 

33 


APPENDIX  B.  OPERATING  INSTRUCTIONS 


A.  PRELIMINARIES 
1.  Hardware  setup 
a.  System  setup 

Setup  system  as  shown  in  Figure  21.  All  signal, 'trigger  connections  should 
be  accompanied  with  a  ground  lead  to  minimize  interference. 


r i bbon 
cable 


Figure  21.  System  interconnection  diagram. 


b.  DASH-16F  switch  setting: 


CHAN  CNFG . 8 

A/D . - . . DIP 

DMA  . 1 


GAIN  (12.5v)  ■ 
Address  (300h) 


dipswitch  A,  B-ON;  C,  D,  USER-OFF 
dipswitch  9,  8--OFF;  7,  6,  5,  4-ON 


c.  Trigger  Requirements 

The  D/A  converter  board  uses  negative  logic  as  shown  in  Figure  22.  The 
pulse  width,  PW,  must  be  less  then  one  pscc  to  ensure  proper  signal  output.  This  is  due 
to  DMA  request  and  acknowledge  hardware.  All  hardware  on  the  correlator  and 
DASH-16F  is  positive  edge-triggered.  The  trigger  magnitude  should  be  approximately 
5.3V  peak,  due  to  CMOS  shift  register  chips  in  the  correlator. 


transai tier 


HH 


PW 


(r~ 

At= 1/6 1440  sec  - ) 

PW<  Ipsec 


+5U 

0 


receiver 
PW  — var  i  ab  1  e 


|€ -  At=  1  /6 1440  sec - »| 

PW 

f - 

Figure  22.  Trigger  specifications. 


d.  Power  Supplies 


+8Vdc  supply . HP6216A  30v  500ma  (or  suitable  substitute) 

— 8Vdc  supply . HP6216A  30v  500ma  (or  suitable  substitute) 

+5Vdc  supply . HP6216A  30v  500ma  (or  suitable  substitute) 


2.  Software  setup 
a.  Transmitter 

1.  Convert  DMAINIT.ASM  and  DMASTOP.ASM  to  BINary  files 
(Ref.  15:  pp.  91-93]., 

2.  Place  files  DMAINIT.BIN,  DMASTOP.BIN,  COMPFFT.1NC,  and 
FFT87B2.INC  in  the  same  directory  as  TRANSMIT.PAS  AND  XMITMES.PAS. 
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3.  Compile  to  disk  TRANSMIT.PAS  and  XMITMES.PAS. 
b.  Receiver 

1.  Compile  TP4D16.PAS  to  disk  to  create  TP4DI6.TPU. 

2.  Place  TP4D16.TPU  in  the  same  director}'  as  RECEIVE.PAS  and  RECMES.PAS. 

3.  Compile  to  disk  RECEIVE.PAS  and  RECMES.PAS. 

B.  MESSAGE  TRANSMISSION  PROCEDURE 

1.  As  necessary  conduct  preliminary  equipment  setup. 

2.  Reset  correlator  (removes  triggers  to  DASH-16F). 

3.  Run  XMITMES  until  "Ready  to  Transmit". 

4.  Run  RECMES  (match  baud  size  used  in  XMITMES). 

5.  When  receiver  prompts  "Ready  to  acquires",  press  enter  key  on  transmitter. 

6.  Demodulated  message  is  displayed  on  receiver  screen  and  stored  in  file 
MESSAGE.DAT. 

TRANSMIT  and  RECEIVE  are  similarly  executed.  However  TRANSMIT  would 
require  manual  loading  of  the  synchronization  baud  in  Table  5  by  individual  loading 
each  tone  with  the  appropriate  symbol. 


Table  5.  SYNCHRONIZATION  BAUD  SYMBOL  SEQUENCE 


tone 

6S 

69 

70 

71 

72 

— 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

symbol 

3 

3 

2 

2 

D 

a 

2 

a 

2 

3 

a 

J) 

D 

a 

3 

1 
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APPENDIX  C.  TESTING  PROCEDURES 


A.  RESPONSE  TESTING  PROCEDURE 

1.  Conduct  preliminary  equipment  setup. 

2.  Run  XMITMES  and  exit. 

3.  Transfer  file  XMITDAT.DAT  to  the  receiver  computer.  (Not  required  for  subse¬ 
quent  response  testing  if  transmitted  message  and  packet  construction  are  un¬ 
changed.) 

4.  Transmit  message  I  AW  message  transmission  procedure.  (RECMES  generates 
output  file  RECSTAT.DAT  used  by  RESPONSE  and  STATS.) 

5.  Run  RESPONSE.  (Ensure  baud  size  is  correctly  edited  into  RESPONSE.) 

6.  Graph  output  file  RESP6au^.DAT. 

B.  SNR  TESTING  PROCEDURE 

1.  Conduct  preliminary  equipment  setup. 

2.  Add  random  noise  source  and  true  RMS  voltmeter  as  shown  in  Figure  23. 

2.2«n 


3.  Set  SNR  level. 

4.  Conduct  steps  1-3  in  the  Response  testing  procedure. 

5.  Run  Statistics  (results  are  in  file  STAT^awfl’.DAT). 
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APPENDIX  E.  TRANSMIT 


program  TRANSMIT; 

(*Transmits  variable  length  packets  up  to  61440.  The  length  is 
determined  by  BCSTARRAY  and  can  be  take  to  a  full  segment  of  memory. 

No  sync  baud  is  used.  All  external  files  are  assumed  to  be  in  the  same 
directory*) 


type 

TNvector  =  array[ 0. . 4095]  of  real; 

TNvectorPtr  =  ATNvector; 

BCSTARRAY  =  array[ -28673. . 32767]  of  byte;  (*stores  output  samples*) 


var 


kx,  (*baud  size*) 

kl,  (*lowest  tone  in  band*) 

k2,  (’"'highest  tone  in  band*) 

NUMBAUDS,  (’’'number  of  bauds  in  packet*) 

BAUDCOUNT,  (-'current  baud  being  processed*) 

BYTECOUNT,  (’"'number  of  bytes  to  be  transferred  by  DMA*) 
I  : integer; 

XREAL,  (*complex  freq/time  domain  arrays*) 

XIMAG  :  TNvectorPtr; 

QAM,  (’-encoding  scheme*) 

INVERSE  :  boolean;  (’'direction  of  FFT*) 

ERROR  :  byte;  (’’'status  of  ComplexFFT*) 

BCST  :  BCSTARRAY;  (’’'packet  storage  buffer*) 

THEFILE  :  file  of  byte;  (’’'time  sequence  output  file*) 

ANSWER  :  char;  (’"'input  variable*) 


(*$I  FFT87B2.  INC*) 

(*$I  COMPFFT.  INC*) 

(* . - . - . - . *) 

procedure  SelectBaud; 

(*SelectBaud  establishes  QPSK  or  QAM  encoding  kx,  kl,  k2,  and  the  number 
of  baud  in  the  packet*) 


var 

ANSWER  : integer; 
MODANSWER  : char; 


begin 

window(20, 10,80,20); 


(•'Select  encoding  scheme*) 
repeat 
clrscr; 

writeln( 1  Select  modulation  desired1); 
writeln; 

write ln( 1  A  —  QPSK'); 

writelnC  B  —  16  QAM'); 

readln(MODANSWER); 
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until  MODANSWER  in  [ ' a'  ' A' , *b' , *B'] ; 
if  MODANSWER  in  ['b\'Bf]  then 
QAM  : =  true; 

(*Select  buad  size*) 
kx:  =0; 
repeat 
c  Irscr* 

if  kx  <  0  then  writeln('TRY  AGAIN1); 
writeln( ' What  is  the  length  of  the  bauds  (kx)?'); 
writeln( 'i. e.  256,  512,  1024,  2048,  4096'); 
readln( ANSWER); 
case  ANSWER  of 
256:  kx:  =256; 

512:  kx:  =512; 

1024:  kx: =1024; 

2048:  kx:  =2048; 

4096:  kx:  =4-  '• 

end;  (*case  AN.  *) 
if  kx  =  0  then  -1; 

until  kx  >  0; 

(*Set  bandlimits*) 
case  kx  of 
256: begin 

kl:  =68;  k2:  =83; 
end; 

512: begin 

kl:  =135;  k2:  =166; 
end; 

1024: begin 

kl:  =269;  k2:=332; 
end; 

2048: begin 

kl:  =537;  k2:=664; 
end; 

4096: begin 

kl:  =1073;  k2: =1328; 
end; 

end;  (*case  kx*) 


(*Select  number  of  baud*) 

write('How  many  bauds  do  your  desire  to  transmit?  '); 
readln(NUMBAUDS); 
window( 1, 1,80,25); 
end;  (*SelectBaud*) 


procedure  DisplayQAM; 

(*DisplayQAM  shows  the  16-QAM  constellation  in  the  upper  portion  of  the 
screen.  A  window  is  set  at  the  bottom  of  the  screen  for  further 


interaction  with  the  user.*) 


begin 

clrscr; 


gotoXY(ll, 1); 
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write(chr(218)); 
for  I:=  12  to  64  do 
write(chr( 196) ); 
write(chr( 191)); 
for  I:=  2  to  19  do 
begin 

GotoXY(ll,I);  write(chr(179)); 

GotoXY(65,I);  write(chr(179)); 
end; 

GotoXY( 11, 20); 
write(chr( 192)); 
for  I  :=  12  to  64  do 
write(chr( 196)); 
write(chr(217)); 


window( 12,2,63,19); 
write('This  program  encodes 


writeln( 

writeln( 

writeln( 


signal. ' ); 

The  vector  space  is 
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!) 
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* 

* 
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* 

* 

) 
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) 
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' ,chr(179) , ' 

) 
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) 
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,chr(179) , 

(*active  window*) 

window( 12,21,80,25); 
end; (*DisplayQAM*) 

(* . *) 

(* . . . *-) 

procedure  DisplayQPSK; 

(*DisplayQPSK  shows  the  QPSK  constellation  in  the  upper  portion  of  the 
screen.  A  window  is  set  at  the  bottom  of  the  screen  for  further 
interaction  with  the  user.*) 


begin 

clrscr; 
gotoXY( 11,1); 
write(chr(218)); 
for  I:=  12  to  64  do 
write(chr( 196)); 
write(chr( 191)); 
for  I:  =  2  to  19  do 
begin 

GotoXY(ll,I);  write(chr( 179)); 

GotoXY(65,I);  write(chr( 179)); 
end; 
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GotoXY( 11,20); 
write(chr( 192) ); 
for  I  :=  12  to  64  do 
write(chr( 196)); 
write(chr(217)); 


window( 12,2,63,19); 


writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 

writeln( 
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writeln( 
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This  program  encodes 
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' ,chr( 179) , 

' ,chr( 179) , 

,chr(179), 

,chr(179), 

,chr( 179) , 

' ,chr( 179) , 1  4 

' ,chr( 179) , ' 


); 

); 

); 

); 

); 

); 

); 

); 

)$ 

); 

); 

); 

); 

); 


(*active  window*) 

window( 12,21,80,25); 
end; (*DisplayQPSK*) 

(* . *) 

(* . *) 

procedure  SelectQAM; 

(*SelectQAM  selects  and  encodes  the  symbols  for  each  tone  in  the  band  set 
by  kl  and  k2.  All  symbols  can  be  randomly  or  manually  selected.*) 


var 

RESPONSE  : char; 

ANSWER, I  : integer; 

VECTORARRAY  :  array[ 0. .  4095]  of  integer; 


(* . *) 

procedure  EncodeData; 

(*EncodeData  loads  the  frequency  domain  arrays  XREAL  and  XIMAG  with  the 
QAM  symbols  in  VECTORARRAY  and  their  complex  conjugate  image.*) 

var 

J  : integer; 

TEMPR, TEMPI  :  real; 

begin 

fillchar(XREALA ,sizeof(XREALA) ,0); 
fillchar(XIMAGA ,sizeof(XIMAGA) ,0); 
for  J:=  kl  to  k2  do 
begin 

case  VECTORARRAY[ J]  of 
0  :  TEMPI;  =  0.0; 

1.  .4  :  TEMPI:  =  90.  0; 

5..  8  :  TEMPI:  =  30.  0; 
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9..  12  :  TEMPI:  =-30.0; 

13. .  16:  TEMPI:  =-90. 0; 
end;-(*case  VECTORARRAY*) 
case  VECTORARRAY[ J]  of 

0  :  TEMPR: =  0. 0; 

1,5,9,13  : TEMPR:  =  90.0 
2,6, 10, 14:  TEMPR:  =  30.0 
3, 7, 11, 15:  TEMPR:  =-30.  0 
4,8, 12, 16:  TEMPR:  =-90. 0 
end; (*case  VECTORARRAY*) 

(*Load  complex  conjugate  image*) 
XREALa[J]  :=  TEMPR; 
XREALa[ kx-J] :  =  TEMPR; 
XIMAGa[J]:=  TEMPI; 
XIMAGa[ kx-JJ :  =  -TEMPI, 
end; (*J:  =kl*) 
end; (*EncodeData*) 


repeat 

clrscr; 

writeln('  Select  one  of  the  following  for  baud  ' ,BAUDCOUNT, ' .  ' ); 
writeln; 

writeln('R  Randomly  select  all  ' ,k2  -  kl  +  1,'  vectors'); 

writeln('l  Individually  select  vectors'); 

readln( RESPONSE); 

until  RESPONSE  in  [ ' r ’ , ' R ' , ' i ? ' I' ] ; 
f il lchar ( VECTORARRAY , s izeof ( VECTORARRAY) , 0) ; 
case  RESPONSE  of 

'r','R'  :  for  I:=  kl  to  k2  do 

VECTORARRAYf I] :  =random( 16)+1; 

* i ' , *  I '  : begin 

I:  =  kl; 

while  I  <=  k2  do 
begin 
repeat 
clrscr; 

writeln( '  Vector  Selection  Menu'); 
writeln('#  vector  desired'); 

writeln( ' 17  random  vector'); 

writeln('l8  no  vector'); 

write( ' 19  no  more  tones  Tone  ',1,'  '); 
readln  (ANSWER); 
until  ANSWER  in  [ 1. . 19] ; 
case  ANSWER  of 

1. .  16: VECTORARRAYf I] :  =  ANSWER; 

17: VECTORARRAYf I] :  =  random( 16)+1; 

18: ; (*nop*) 

19: 1:  =  k2; 
end; (*case  ANSWER*) 

I:  =1+1; 

end; (*while  I*) 
end; (*' i ' , ' I '  *) 
end; (*case  RESPONSE*) 


EncodeData; 
end; (*SelectQAM*) 


(* . *) 

procedure  SelectQPSK; 

(*SelectQPSK  selects  and  encodes  the  symbols  for  each  tone  in  the  band 
set  by  kl  and  k2.  All  symbols  can  be  randomly  or  manually  selected.*) 

var 

RESPONSE  : char; 

ANSWER, I  : integer; 

PHASEARRAY  : array[ 0. . 4095]  of  integer; 

(* . *) 

procedure  EncodeData; 

(*EncodeData  loads  the  frequency  domain  arrays  XREAL  and  XIMAG  with  the 
QPSK  symbols  in  VECTORARRAY  and  their  complex  conjugate  image. *) 


J  : integer; 

TEMPR, TEMPI  .-real; 

begin 

fillchar(XREALA,sizeof(XREALA) ,0); 
f illchar(XIMAGA ,sizeof(XIMAGA) ,0); 
for  J:=  kl  to  k2  do 
begin 

case  PHASEARRAY[  J]  of 
0  :  TEMPI:  =  0.0; 

1.2  :  TEMPI:  =  80.0; 

3.4  :  TEMPI:  =-80.0; 
end; (*case  PHASEARRAY*) 
case  PHASEARRAY[ J]  of 

0  :  TEMPR:  =  0.  0; 

1.4  :  TEMPR:  =  80.0; 

2.3  : TEMPR:  =-80.  0; 
end; (*case  PHASEARRAY*) 
XREALa] J]  :  =  TEMPR; 

XREALa] kx-J)  :=  TEMPR; 

XIMAG  a  (  J):=  TEMPI; 

XIMAGa[ kx-J] :  =  -TEMPI; 

end; (*J:  =kl*) 
end; (*EncodeData*) 


begin 

repeat 

clrscr; 

writeln( '  Select  one  of  the  following  for  baud  ' , BAUDCOUNT , 1 .  1 ) ; 
writeln; 

writeln('R  Randomly  select  all  1 ,k2  -  kl  +  1,'  phases'); 
writeln('l  Individually  select  phases'); 

read ln( RESPONSE); 

until  RESPONSE  in  [ ' r' , 'R' , ' i' , ' I ' ] ; 
f i 1 1 cha  r ( PHASEARRAY , s izcof ( PHASE ARRAY ) , 0 ) ; 
case  RESPONSE  of 

' r ' , ' R' : for  I:  =  kl  to  k2  do 
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PHASEARRAY[ I] :  =random(4)+l; 

' i' , *  I  * :  begin 

I:  =  kl; 

while  I  <=  k2  do 
begin 
repeat 
clrscr; 

writeln( 1  Phase  Selection  Menu'); 
writeln('#  phase  desired'); 

writeln('5  random  phase' ); 

writeln('6  no  phase'); 

write('7  no  more  tones  Tone  ',1,'  '); 
readln  (ANSWER); 
until  ANSWER  in  [ 1. . 7] ; 
case  ANSWER  of 

1. . 4: PHASEARRAY[ I] :  =  ANSWER; 

5: PHASEARRAY[ I] : =  random(4)+l; 

6: ; (*nop*) 

7:  I:  =  k2; 

end; (*case  ANSWER*) 

I:  =1+1; 

end;(*while  I*) 
end; (*' i' , ' I'*) 
end; (*case  RESPONSE*) 

EncodeData; 
end; (*SelectQPSK*) 


procedure  ScaleData; 

(*ScaleData  converts  each  real  value  in  array  XREAL  down  to  a  byte  and 
stores  the  byte  in  packet  storage  buffer  BCST.  INDEX  sets  the  location 
in  the  buffer  of  each  byte*) 


INDEX, J, TEMP  : integer; 

DATA  : byte; 


begin 

for  J  :  =  0  to  kx-1  do 
begin 

TEMP  :=  round(XREALA[ J]  +  126); 
if  TEMP  <  0  then 
TEMP  :=  0; 

DATA  :  =  TEMP’ 

INDEX  :=  J  +’ (BAUDCOUNT  -  1)  *  kx  -  28673; 

BCST[ INDEX]  :  =  DATA; 
write( THEFILE , DATA) ; 
end;  (*for  J*) 
end;  (*ScaleData*) 

(* . *) 

(•" . - . *) 

procedure  Dmainit(var  BCST  :  BCSTARRAY; 

BYTECOUNT  : integer); 

(''•'Assembly  language  procedure  used  to  initialize  and  unmask  the  DMA  for 
data  transfer.  The  source  must  be  converted  to  a  BIN  file.*) 
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external  'DMAINIT5.BIN'; 

(* . *) 

procedure  Dmastop;  (*Masks  DMA,  stopping  data  transfer*) 
external  ’DMASTOP.BIN'; 

(* . *) 

procedure  TINPUT; 

(*This  procedure  reads  in  binary  file,  H0UTFA1.DAT,  into  an  array, 
performs  an  FFT  using  ComplexFFT  and  displays  the  results 
with  GraphData*) 

var 

NUMPTS  :  integer; 

(* . *) 

procedure  GetData; 

(*GetData  reads  time  domain  sequence  from  file  H0UTFA1.DAT  into  XREAL*) 
var 

THEFILE  : file  of  byte; 

DATA  : byte; 

I, 

TEMP  : integer; 

rdata  : real; 

begin 

as sign( THEFILE,  'H0UTFA1.  DAT' ); 

reset(THEFILE); 

new( XREAL); 

new(XIMAG); 

fil lchar( XREALa ,sizeof  (XREALa) ,0); 
f illchar(XIMAGA ,sizeof  (XIMAGa) ,0); 
for  I  : =  0  to  kx-1  do 
begin 

read( THEFILE, DATA);  (*read  file  one  byte  at  a  time*) 

TEMP  :=  DATA; 

temp  :  =  temp  and  255;  (*reduces  number  of  bits  used*) 

(*to  represent  the  input  data*) 

XREALa[I]  :=  (TEMP  -  126);  (*loads  XREAL  array*) 
end;  (*for  I*) 
close( THEFILE); 
end;  (*GetData*) 


(* . - . - . *) 

(* . *) 

procedure  PutData; 

(*Writes  frequency  domain  arrays  XREAL  and  XIMAG  to  file  TFFTOUT.  OUT*) 
var 

I  : integer; 


OUTFILE  : text; 
begin 

ass ign( OUTFILE , ' TFFTOUT. DAT ’ ) ; 
r  ewr it  e ( OUTF I LE ) ; 
for  I  :  =  0  to  kx-1  do 

writeln( OUTFILE, I, ’  ' , XREALa [ I] , 1  ’ ,XIMAGa[ I] ); 
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close(OUTFILE); 
end;  (*PutData*) 

(* . *) 

procedure  GraphData; 

(*GraphData  plots  the  kx/2  points  of  complex  array  XREAL  and  XIMAG. 
Symbols  are  decoded  and  represented  by  the  color  of  each  tones 
spectral  line.*) 

var 

I,J,  (*loop  indices*) 

KXDIV2, 

COLOR, 

X, XZ, 

Y, YZ  : integer; 

begin 

GraphColorMode; 

palette(2); 

Gr aphbackground ( 0 ) ; 

draw(30, 10,30, 170, 1);  (*y-axis*) 

draw(30, 170,300, 170, 1);  (*x-axis  for  full  spectrum*) 

draw( 30, 155, 300, 155,1); 

draw(30,80,300,80,l);  (*x-axis  for  zoom  spectrum*) 

(*every  other  dot  equals  one  bin  at  256*) 
draw( 166, 170, 166, 174,1);  (*x-axis  scale  marks*) 

draw( 196, 170, 196, 174, 1);  (*x-axis  scale  marks*) 

draw(30, 113,26, 113,1);  (*y-axis  scale  marks*) 

draw(30,23,26,23,l);  (*y-axis  scale  marks*) 

I:  =  23; 

(*Grouping  of  4  horizontal  lines  to  identify  background  color 
used  for  one  of  the  symbols*) 
repeat 

draw(31, I, 300,1,0); 
draw(31, 1+1,300, 1+1,1); 
draw(31,I+2,300,I+2,2); 
draw( 31, 1+3, 300,1+3, 3); 

I:  =1+12; 
until  I  >  75; 

KXDIV2  :=  round(kx  /  2); 


(*Decode  svmbols  and  assign  color*) 
for  I  :=  0  TO  KXDIV2  do 
begin 

if  (XREALa[I]  >=  0)  and  (XIMAGa[I)  >  0)  then 
COLOR  :=  0; 

if  (XREALa[I]  <  0)  and  (XIMAGa[IJ  >  0)  then 
COLOR  :=  1; 

if  (XREALa[I]  <  O')  and  (XIMAGa[I]  <=  0)  then 
COLOR  :=  2; 


if  (XREALa[I]  >=  0)  and  (XIMAGa[I]  <=  0)  then 


COLOR  :=  3 


(*zoom  spectrum*) 


if  (I  >=  kl)  and  (I  <=  k2)  then 
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begin 

XZ  :=  round(((I  *  4096.0)  /  kx)  -  1040); 

YZ  :  =  80  -  round(0.  5*(sqrt(sqr(XREALA[ I]  )  +  sqr(XIMAG/.[  I]  )))); 
dr aw( XZ , 80 , X? , YZ , COLOR) ; 
end; 

(*full  spectrum*) 

X  :=  round( I/KXDIV2  *  256  +30); 

Y  : =170-round(20*ln(sqrt(sqr(XREALA[  I]  +1)+ 
sqr(XIMAGA[I]+l)))/ln(10)); 
draw(X, 170, X,Y, COLOR); 
end; 

end;  (*GraphData*) 

(* . *) 

begin 

INVERSE  :=  false;  (*sets  forward  FFT*) 

ERROR  :=  0; 

GetData; 

NUMPTS  :=  kx; 

writeln( 1  Performing  FFT1); 

Comp lexFFT  ( NUMPTS , INVERSE , XREAL , XI MAG , ERROR ) ; 
writeln( 'Error  =  1 ,ERR0R,'  hit  the  enter  key1 );  readln; 

PutData; 


GraphData; 
end;  (*Tinput*) 

(* . . *) 

(* . *) 


begin 

repeat 

clrscr; 

assign(THEFILE, 'H0UTFA1. DAT' ); 

rewrite(THEFILE); 

new( XREAL); 

new(XIMAG); 

INVERSE: =true; 

QAM:  =false; 

SelectBaud; 

if  QAM  then 
DisplayQAM 
else 

DisplayQPSK; 

(*Packet  contruction  loop*) 

for  baudcount  :  =  1  to  numbauds  do 
begin 

if  QAM  then 
SelectQAM 
else 

SelectQPSK; 

writeln( ' Performing  IFFT’); 
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ComplexFFT(kx, INVERSE ,XREAL , XIMAG , ERROR) ; 
ScaleData; 

end;  (*for  BAUD COUNT*) 

dispose(XREAL); 

dispose(XIMAG); 

close(THEFILE); 

BYTECOUNT  :=  NUMBAUDS*kx-l; 

writeln( 'Press  return  to  transmit' ); readln; 

Dmainit(BCST, BYTECOUNT); 

repeat 

writeln( 'Transmit  some  more?  (*yes  or  no*)  '); 
readln( ANSWER); 

until  ANSWER  in  [ V , 'N' , *y‘ , ' Y* ] ; 
dmastop; 

window( 1, 1,80,25); 
until  ANSWER  in  [  'n' ,'N']; 

Tinput; 

end. 
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APPENDIX  F.  XMITMES 


program  XMITMES; 

( ‘'Transmits  a  syncbaud  and  message  from  file  'MESSAGE.DAT1.  The  message 
is  encoded  using  QPSK.  'MESSAGE.DAT1  is  a  text  file.  It  should  already 
exist  before  using  this  program.  Output  is  used  to  collect  data  for 
TESTING*) 

const 

FIRST_ELEMENT  =  -28929; 
type 

TNvector  =  arrayf  0. . 4095]  of  real; 

TNvectorPtr  =  ATNvector; 

(*Sync  +  61440*) 

BCSTARRAY  =  arrayf FIRST_ELEMENT. .  32767]  of  byte; 


var 


kx, 

kl,k2, I ,w, 

NUMBAUDS , MAXNUMBAUDS , 
BAUDCOUNT, BYTECOUNT, 
SYMBOLCOUNT , MAXNUMCHAR , 


MESSAGES I 2E 
MAGNITUDE, 
CHARACTERS_PER_BAUD 
XREAL,XIMAG 
INVERSE 

TEMPBYTE, ERROR 

BUST 

BYTEFILE 

TESTFILE 

ANSWER, 

NEXTCHAR 

(*$I  FFT87B2.  INC*) 

(*$I  COMPFFT. INC*) 


: integer; 

:  real; 

: TNvectorPtr; 

: boolean; 

: byte; 

: BCSTARRAY; 

: file  of  byte; 
: text; 

: char; 


(* . *) 

procedure  SyncBaud; 

(••''Process  the  synchronization  baud  and  stores  the  256  point  time  domain 
sequence  at  the  beginning  of  the  packet  storage  area.*) 


var 

SYNCREAL,  SYNCIMAG  :  TNvectorPtr; 

J,  TEMP  :  integer; 

SYNCDATA  :  byte; 

SYNCMAG  :  real; 

begin 

new( SYNCREAL); 
new( SYNCIMAG); 

SYNCMAG:  =  MAGNITUDE; 
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(*load  synchronization  symbols*) 
f i 1 1 char ( SYNCREALa , s izeof ( SYNCREALa ) , 0 ) j 
f illchar( SYNCIMAGa , sizeof ( SYNCIMAGa ) , 0) ; 

SYNCREALa [68]:=  -SYNCMAG;  SYNCIMAGa [ 68]  :  =  -SYNCMAG 
SYNCREALa [69]:=  -SYNCMAG;  SYNCIMAGa [ 69] :  =  -SYNCMAG 
SYNCREALa [ 70] : =  -SYNCMAG;  SYNCIMAGa[  70]  :  =  SYNCMAG 
SYNCREALa [ 71]:=  -SYNCMAG;  SYNCIMAGa [  71]  :  =  SYNCMAG 
SYNCREALa [72]:=  SYNCMAG;  SYNCIMAGaJ 72] :  =  -SYNCMAG 
SYNCREALa [  73] :  =  SYNCMAG;  SYNCIMAGa[  73] :  =  SYNCMAG 
SYNCREALa [74]:=  -SYNCMAG;  SYNCIMAGa[ 74] : =  SYNCMAG 
SYNCREALa [ 75]  :  =  SYNCMAG;  SYNCIMAGa [ 75] : =  SYNCMAG 
SYNCREALa [ 76] : =  -SYNCMAG;  SYNCIMAGa [  76]  :  =  SYNCMAG 
SYNCREALa[ 77]:=  -SYNCMAG;  SYNCIMAGa [ 77] : =  -SYNCMAG 
SYNCREALa [ 78]  :  =  SYNCMAG;  SYNCIMAGa[ 78] : =  -SYNCMAG 
SYNCREALa[  79] :  =  -SYNCMAG;  SYNCIMAGa[ 79] : =  -SYNCMAG 
SYNCREALa [  80]  :  =  SYNCMAG;  SYNCIMAGa[ 80]  :  =  SYNCMAG 
SYNCREALa [  81] :  =  SYNCMAG;  SYNCIMAGa[ 81] :  =  -SYNCMAG 
SYNCREALa [ 82]  :  =  -SYNCMAG;  SYNCIMAGa [ 82] : =  -SYNCMAG 
SYNCREALa [83]:=  SYNCMAG;  SYNCIMAGa[ 83] : =  SYNCMAG; 

(*complex  conjugate  image*) 
for  J  :=  68  to  83  do 
begin 

SYNCREALa [ 25 6 -J]  :  =  SYNCREALa] J] ; 

SYNCIMAGa [ 256-J] : =-SYNCIMAGa[ J] ; 
end;  (*for  J*) 

ComplexFFT( 256 , INVERSE , SYNCREAL, SYNC IMAG .ERROR) ; 

(*scale  data/ load  time  sequence*) 
for  J  : =  0  to  255  do 
begin 

TEMP: =round( SYNCREALa [J]  +  126); 

if  TEMP  <  0  then 
TEMP:  =0; 

SYNCDATA:  =TEMP; 

BCST[  J+FIRST_ELEMENT]  :  =SYNCDATA; 
end; (*for  J*) 
disposed SYNCREAL); 
dispose( SYNCIMAG) ; 
end; (*SyncBaud*) 


procedure  SelectBaud; 

(*SelectBaud  establishes  kx,  kl,  and  k2*) 
var 

ANSWER  : integer; 

begin 
kx:  =0; 

(*select  baud  size*) 
repeat 

if  kx  <  0  then  writeln('TRY  AGAIN'); 

writelnC 1  What  is  the  length  of  the  bauds  (kx)?'); 
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writeln( ' i. e.  256,  512,  1024,  2048,  4096'); 
readln( ANSWER); 
case  ANSWER  of 

256:  kx:  =256; 

512:  kx:=512; 

1024:  kx: =1024; 

2048:  kx: =2048; 

4096:  kx: =4096; 
end;  (*case  ANSWER*) 
if  kx  =  0  then  kx  :=  -1; 
until  kx  >  0; 


(*set  tone  limits*) 
case  kx  of 
256: begin 

kl:  =68;  k2:=83; 
end; 

512:  begin 

kl:  =135;  k2:  =166; 
end; 

1024: begin 

kl:  =269;  k2:  =332; 
end; 

2048: begin 

kl:  =537;  k2:=664; 
end; 

4096: begin 

kl:  =1073;  k2: =1328; 
end; 

end;  (*case  kx*) 
end;  (*SelectBaud*) 

(* . - . *) 

(* . - . *) 

procedure  TailorPacket; 

(*TailorPacket  sets  the  maximum  number  of  baud  required  to  encode  the 
message*) 


begin 

MESSAGESIZE:  =  filesize(BYTEFILE); 

writeln( 'Message  is  1 , MESSAGESIZE, '  bytes.'); 

(*kx/2  is  the  number  of  bit/baud  for  QPSK.  kx/2-2  is  the  number 
for  DQPSK.  Each  character  is  8  bits*) 

CHARACTERS_PER_BAUD: =(kx/8  -  2)/8; 

(*61440/kx  is  maximum  number  of  bauds  possible*) 

MAXNUMCHAR:  =  trunc(61440.  0/kx  *  CHARACTERS_PER_BAUD); 
if  MESSAGESIZE  >  MAXNUMCHAR  then 
begin 

writeln( 'Message  is  to  large.  The  last  ', 

MESSAGESIZE  -  MAXNUMCHAR, 

'  characters  will  not  be  transmitted.'); 
MESSAGESIZE: =MAXNUMCHAR; 
end; 

MAXNUMBAUDS: =trunc(MESSAGESIZE  /  CHARACTERS_PER_BAUD); 
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(*ensure  last  few  characters  are  included*) 
if  frac(MESSAGESIZE  /  CHARACTERS  PER_BAUD)  >0.0  then 
MAXNUMBAUDS: =MAXNUMBAUDS  +  1; 

repeat 

writeln; 

writeln( 'Enter  number  of  ',kx,'  bauds  to  process.  ', 

MAXNUMBAUDS , '  is  the  maximum.'); 
readln(NUMBAUDS); 

until  NUMBAUDS  in  [ 1. .  MAXNUMBAUDS) ; 
end;  (*TailorPacket*) 

(* . *) 

(* . - . *) 

procedure  DiffEncode; 

(*DiffEncode  differential  encodes  symbols  on  a  tone-to-tone  basis. 
BYTEFILE  is  read  from  one  byte  at  a  time.  The  byte  is  isolated  into 
2-bit  groups  and  stored  in  BITS.  BITS  is  then  used  to  DQPSK  encode  the 
frequency  domain  arrays  XREAL  and  XIMAG.  Bytes  partially  encoded  are 
carried  over  into  the  next  baud  by  global  variable  TEMPBYTE.*) 

var 

J  : integer; 

BITS  :byte; 

begin 

fillchar(XREALA,sizeof(XREALA),0); 

fillchar(XIMAGA,sizeof(XIMAGA),0); 

(*first  tone  of  every  baud  set  to  pi/2*) 

XREALa[ kl] :  =  MAGNITUDE; 

XIMAGa[ kl]  :  =  MAGNITUDE; 

if  SYMBOLCOUNT  =  0  then 

read( BYTEFILE , TEMPBYTE); 

(*break  apart  character*) 
for  J:=  (kl  +  1)  to  k2  do 
begin 

SYMBOLCOUNT:  =SYMB0LC0UNT  +  1; 
if  frac( SYMBOLCOUNT  /  4)  =  0.25  then 
BITS: =  (TEMPBYTE  and  $C0)  shr  6; 
if  frac( SYMBOLCOUNT  /  4)  =  0. 5  then 
BITS: =  (TEMPBYTE  and  $30)  shr  4; 
if  frac( SYMBOLCOUNT  /  4)  =  0.  75  then 
BITS:  =  (TEMPBYTE  and  $0C)  shr  2; 
if  frac( SYMBOLCOUNT  /  4)  =  0. 0  then 
begin 

BITS: =  TEMPBYTE  and  $03; 
if  not  E0F( BYTEFILE)  then 
read( BYTEFILE, TEMPBYTE) 
else 

TEMPBYTE:  =$40;  (*f ill  character*) 

if  (BITS  <  0)  and  (BITS  >  3)  then 

writeln( ' Bits  not  assigned  properly'); 


(■^differential  encode*) 
case  BITS  of 

0:  begin  XREALa[J]:=  XREALa[J-1]; 

XIMAGa[ J] :  =  XIMAGa[ J-l] ;  end; 
1:  begin  XREALa[ J] :  =-XIMAGa[ J-l] ; 

XIMAGa [  J] :  =  XREALa [ J-l] ;  end; 
2:  begin  XREALa[J]:  =  XIMAGa  [J-l]; 

XIMAGa [  J] : =-XREALa[  J-l] ;  end; 
3:  begin  XREALa[ J] j =-XREALa[ J-l] ; 

XIMAGa[ J] : =-XIMAGa[ J-l] ;  end; 
end; (*case  BITS*) 
end;  (*for  J*) 


(*complex  conjugate  image*) 
for  J:=  kl  to  k2  do 
begin 

XREALa  [  kx  -  J]:=  XREALa  [J]; 

XIMAGa [ kx  -  J]:=-XIMAGa[J]; 

writeln(TESTFILE , BAUDCOUNT: 4 , J: 5 , trunc( XREALa [  J]  ) :  6 , 
trunc( XIMAGa [ J] ) : 6 ) ; 

end; 

end; (*Diff Encode*) 

. - . *) 

(* . - . . . . . *) 

procedure  ScaleData; 

(*ScaleData  converts  each  real  value  in  XREAL  down  to  a  byte  and  stores 
the  byte  in  the  packet  storage  buffer,  BCST.  INDEX  establishes  the 
location  in  the  buffer  of  each  byte  in  the  packet.*) 


var 

INDEX, J, TEMP  : integer; 

DATA  : byte; 

begin 

for  J  :  =  0  to  kx-1  do 
begin 

TEMP  :=  round( XREALa [J]  +  126); 
if  TEMP  <  0  then 
TEMP  :=  0; 

DATA  :=  TEMP; 

(*256  is  added  to  INDEX  to  start  message  bauds 
after  the  sync  baud*) 

INDEX  :=  J+( BAUDCOUNT- l)*kx+FIRST_ELEMENT+25 6; 

BCST[ INDEX]  :  =  DATA; 
end; (*for  J*) 
end;  (*ScaleData*) 

(*--- - - - - - - - *) 

(* . *) 

procedure  Dmainit(var  BCST  rBCSTARRAY; 

BYTECOUNT  : integer); 

(*Assembly  language  procedure  used  to  initialize  and  unmask  the  DMA  for 
data  transfer.  The  source  code  must  be  converted  to  a  BIN  file.*) 


external  'DMAINIT5.BIN'; 

(*— - . . - . *) 
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procedure  Dmastop; 

(*Masks  DMA,  stopping  data  transfer.*) 
external  'DMASTOP.  BIN1; 

C* . *) 

(* . *) 

begin 

clrscr; 

(*contains  hex  values  to  be  encoded  and  transmitted*) 
assign(BYTEFILE, 'MESSAGE.  DAT' ); 
reset(BYTEFILE); 

(’■'Output  file  of  encoded  symbols.  Used  for  system  testing*) 
assign(TESTFILE, 'XMITDAT.  DAT’ ); 
rewrite(TESTFILE); 

INVERSE:  =true; 
repeat 

writeln( 'Enter  magnitude  of  tones,  (must  be  >  0.0)  '); 
read ln( MAGNITUDE); 
until  MAGNITUDE  >  0.  0; 

writeln( 'Loading  sync  baud. '); 

SyncBaud; 

SelectBaud; 

TailorPacket; 

SYMBOLCOUNT:  =0; 

TEMPBYTE: =$00; 

writeln( ' Number  of  bauds  is  ' ,numbauds); 

new(XREAL); 

new(XIMAG); 

for  baudcount  : =  1  to  numbauds  do 
begin 

DiffEncode; 

writeln( ' Performing  IFFT  ' , BAUDCOUNT, '  ' , 

NUMBAUDS -BAUDCOUNT,'  left’); 

Comp 1 exFFT( kx , INVERSE , XREAL , XIMAG , ERROR) ; 

ScaleData; 

end;  (*for  BAUDCOUNT’’') 

dispose(XREAL); 

dispose(XIMAG); 

close(BYTEFILE); 

close(TESTFILE); 

BYTECOUNT  :=  256  +  NUMBAUDS*kx  -  1; 
repeat 

writeln( 1  Press  return  to  transmit’ );  readln; 

Dmainit( BCST, BYTECOUNT) ; 
repeat 

writeln( 'Transmit  some  more?  (*ves  or  no*)  '); 
readln(ANSWER); 

until  ANSWER  in  [ 'n' , ' N' , ' y* , ’ Y' ]  ; 
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dmastop; 

until  ANSWER  in  [  V/N']; 

(*  reset(TESTFILE); 

while  not  EOF(TESTFILE)  do 
begin 

while  not  EOLN(TESTFILE)  do 
begin 

read(TESTFILE,NEXTCHAR); 

write(NEXTCHAR); 

end; 

readln(TESTFILE); 

writeln; 

end; 

close(TESTFILE);  *) 
end. 
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APPENDIX  G.  DMAINIT 


codeseg  segment 

public  dmainit 
assume  cs: codeseg 

j procedure  DMAINIT  (  BLKADDRESS  :  XMITPOINTER; 

;  BYTECOUNT  :  INTEGER); 

; this  procedure  initializes  dma  channel  3  and  sets  the 
; parameters  to  output  the  array  best  by  passing  the  address 
; of  the  array  start  on  the  stack.  BYTECOUNT  is  the  number 
; of  bytes  to  transfer  and  is  pushed  on  the  stack  by  the 
; calling  program 

dma  equ  0 

dmapage  equ  80h 

dmainit  proc  near 


push 

bp 

mov 

bp.sp 

;use  bp  to  address  stack 

les 

di.dword  ptr[bp+6j 

;move  add  of  best  into  es:  di 

mov 

al,5bh 

; dma  chan  3  single  mode,  read,  autoinit 

out 

dma+11, al 

out 

dma+12,al 

; reset  first/last  ff 

mov 

ax,es 

; calc  high  order  4  bits  of  buffer  area 

mov 

cl, 4 

rol 

ax, cl 

push 

ax 

; save  ax  for  dma  start  addr 

and 

al,0fh 

out 

dmapage+2,al 

; store  in  ch  3  dma  page  reg 

pop 

ax 

and 

al,0f0h 

add 

ax,di 

; get  page  offset 

out 

dma+6,al 

; output  waveform  buffer  start  addr 

mov 

al,ah 

out 

dma+6,al 

mov 

ax,[ bp+4] 

; output  dma  byte  count 

out 

dma+7 ,al 

mov 

al  ,ah 

out 

dma+7 , al 

mov 

al,3 

; unmask  ch  3  to  start 

out 

dma+10,al 

pop 

bp 

ret 

6 

; pop  6  bytes  off  stack  for  aadr  of  best 

dmainit  endp 
codeseg  ends 


APPENDIX  H.  DMASTOP 


codes eg  segment 

public  dmastop 
assume  cs: codes eg 

; procedure  DMASTOP; 

) 

; this  procedure  stops  dma  channel  3 

dma  equ  0 

dmapage  equ  8 Oh 

dmastop  proc  near 

push  bp 

mov  al,7  ;mask  ch  3  to  stop 

out  dma+10,al 

pop  bp 

ret 

dmastop  endp 
codes eg  ends 
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APPENDIXI.  RECEIVE 


program  RECEIVE; 

(*This  program  acquires  an  analog  signal,  stores  the  raw  data  in  a  memory 
buffer,  converts  the  data  to  TYPE  real,  performs  an  FFT  using  COMPFFT, 
differential  decodes  between  adjacent  tones  and  displays  the  result  in 
the  frequency  domain  with  the  phase  quadrant  represented  by  color. 

Data  input  can  also  be  read  in  from  the  time  domain  sequence  file 
H0UTFA1.DAT  generated  by  program  TRANSMIT*) 

uses  Graph,  Crt,  tp4dl6; 

(*$I-*) 

(*$R-*) 


const 

kx  =  256;  (*set  baudlength*) 
kxml  =255; 

type 

TNvector  =  array[ 0. .  kxml] 
TNvectorPtr  =  ATNvector; 


var 

INVERSE  : boolean; 

XREAL,  XIMAG  : TNvectorPtr; 

ERROR  : byte; 

NUMPTS, 

kl,k2  : integer; 

ANSWER  : char; 

(*$I  FFT87B2.  INC*) 

(*$I  COMPFFT.  INC*) 

(*- . *) 

(* . *) 


procedure  AcquireData; 

(*AcquireData  intializes  Metrobyte  DASH  -16F  data  acquisition  board. 
Using  TTOQLS  procedures  D16_init  and  D16_ainm.  Data  transfer  is 
controlled  by  the  DMA  controller  and  initialized  by  D16_ainm,  and 
disabled  by  D16_dma_int_disable.  TT00LS  procedures  are  external 
procedures  included  by  'uses'  tp4dl6.*) 

const  max_buffer  =  1000; 

var  i:  integer; 

rate:  real; 

cnt_num,  mode,  cycle,  trigger, 
base_adr,  err_code,  int_level,  droa_level, 
board_num,  chanlo, 

op_type,  status,  next_cnt,  err_code_s  :  integer; 


of  real; (*TYPE  for  real  and  imaginary 
data  for  FFT  routing*) 

(•'Pointer  for  FFT  data  array 
which  allows  dynamic  allocation 
of  memory*) 
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dmaPointer:  pointer; 
datavector:  a  integer; 

ad_data,  chan_data:  array[0. . max_buffer]  of  integer; 

begin 

clrscr; 

(’’'allocates  memory  on  the  heap  and  returns  a  pointer  to  the  start 
of  the  buffer*) 

GetDMABuf  f er (max_buf f er , dmaPointer , err_code) ; 

(*This  statement  assigns  a  generic  pointer  to  a  variable  of  a 
specific  pointer  type,  i. e.  Ainteger,  so  that  the  pointer  can  be 
passed  to  the  dl6__ainm  routine.*) 
datavector  :  =  dmaPointer; 

board_num  :=  0;  int_level  :=  7;  dma_level  :=  1; 
base_adr  :  =  $300; 

(*initialize  the  driver*) 

dl6_init ( board_num , base_adr , int_level , dma_level , err_code) ; 
chanlo  :  =  0; 

cycle:  =0;  (*0  -  one  sweep  of  the  DMA  1  -  autoinitialize*) 

trigger: =0;  (*0  -  external  1  -  internal*) 

cnt_num:  =kx;  (*#  of  samples*) 

rate  :  =  10000.0;  (*used  for  internal  trigger*) 

mode  :=  2;  (*DMA  mode*) 

writeln( 'Ready  to  acquire'); 

(*colects  kx  analog  values  using  DMA  and  stores  in  a  buffer*) 
dl6_ainm( board_num , chanlo , mode , cycle , trigger , cnt_num , rate , 
datavector a , err_code) ; 

(*status  indicates  the  progress  of  acquisition.  When  all 
samples  have  been  acquired  status=0*) 
status  :=  11; 

(*wait  until  all  data  acquired*) 
repeat 

dl6_dma_int_status ( board_num , op_type , status , next_cnt ,  err_code_s ) ; 
until  status  =  0; 
writeln( 'Data  received'); 

if  err_code  <>  0  then 

dl6_print_error( err_code) 
else 
begin 

writelnC ' Processing  data'); 

(*converts  left  justified  data,  returns  the  true  binary  value 
of  the  sampled  data*) 

dl6_convert_data(2047 ,cnt_num, datavector a ,ad_data[ 0] , 
chan_data[ 0] ,0,err_code); 

new(XREAL); 

new(XIMAG); 

fillchar(XREALA,sizeof  (XREALa) ,0); 
fillchar(XIMAGA,sizeof  (XIMAGa) ,0); 
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for  I:=  0  to  kxml  do 
begin 

XREALa[I]  :=  ad_data[  I]  /5; 
end; 

end; 

(*stop  DMA  operation  if  in  autoinitialization*) 
dl6_dma_int_disable(  board_nutn ,  er r_code) ; . 

(*frees  memory  allocated  with  GetDMABnffer*) 

FreeDMABuf  f er( max_buf f er , dmaPointer , err_code) ; 

end; (*Acquire*) 

(* . *) 

c* . - . *) 

procedure  GetData; 

(*GetData  loads  XREAL  with  the  real  time  domain  sequence  generated  by 
program  TRANSMIT.  Only  one  baud  is  loaded*) 

var 

THEFILE  : file  of  byte; 

DATA  : byte; 

I, 

TEMP  : integer; 


begin 

assign(THEFILE, 'H0UTFA1.  DAT' ); 

reset(THEFILE); 

new( XREAL); 

new(XIMAG); 

fillchar(XREALA,sizeof  (XREALa) ,0); 
fillchar(XIMAGA,sizeof  (XIMAGa),0); 
for  I  :=  0  to  kxml  do 
begin 

read(THEFILEjDATA);  (*read  file  one  byte  at  a  time*) 

TEMP  :=  DATA;  (*puts  byte  into  integer  variable*) 

XREALa[I]  :=  (TEMP  -  126);  (*loads  XREAL  array*) 
end;  (*for  I*) 
close(THEFILE); 


end;  (*GetData*) 

(* . 

(* . 


*) 

•*) 


procedure  DiffDecode; 

(*DiffDecode  differentially  decodes  complex  frequency  domain  arrays  XREAL 
and  XIMAG.  Four  decoded  symbols  are  recombined  into  a  byte  and 
transferred  to  file  BYTES0UT.DAT.*) 


var 

I, 

SYMBOLCOUNT  : integer; 

TEMPREAL , TEMPIMAG  : real; 
BITS,TEMPBYTE  : byte; 

TEMPCHAR  : char; 

OUTFILE  :  TEXT; 

begin 
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as sign( OUTFILE, 'BYTESOUT.  DAT'); 
rewr ite( OUTFILE ) ; 

SYMBOLCOUNT:  =0; 

TEMPBYTE:  =0; 

for  I  :=  kl  to  (k2  -  1)  do 
begin 

(*Complex  multiply  two  adjacent  tones,  I  and  the  complex  conjugate  of 
1+1.  This  will  give  the  phase  difference  between  the  two  tones. 

The  answer  is  in  rectangular  notation*) 

TEMPREAL:  =XREALa[ I]  *  XREALa[  1+1]  +  XIMAGa[I]  *  XIMAGa[ 1+1] ; 
TEMPIMAG: =XREALa[ I]  *  XIMAGa[ 1+1]  -  XREALa[ 1+1]  *  XIMAGa[ I] ; 

(*Complex  multiply  (TEMPREAL  +  j  TEMPIMAGE)  and  (1+j).  This  rotates 
the  differential  vector  45  degrees.  XREAL  and  XIMAG  are  used  to 
store  the  results.  This  eliminate  the  original  data*) 

XREALa[ I] :=( TEMPREAL  -  TEMPIMAG)/80; 

XIMAGa[ I] :=( TEMPREAL  +  TEMPIMAG)/80; 

if  (XREALa[I]  >=  0)  and  (XIMAGa[ I]  >  0)  then 
BITS:  =$00; 

if  (XREALa[I]  <  0)  and  (XIMAGa[I]  >  0)  then 
BITS:  =$01; 

if  (XREALa[I]  <  0)  and  (XIMAGa[ I]  <=  0)  then 
BITS:  =$03; 

if  (XREALa[I]  >=  0)  and  (XIMAGa[I)  <=  0)  then 
BITS: =$02; 

SYMBOLCOUNT  :  =  SYMBOLCOUNT  +  1; 

(*f ill  TEMPBYTE  with  four  symbols*) 
if  frac( SYMBOLCOUNT  /  4)  =  0.25  then 
TEMPBYTE  :=  (BITS  shl  6); 
if  frac( SYMBOLCOUNT  /  4)  =  0.  5  then 
TEMPBYTE  :=  (BITS  shl  4)  or  TEMPBYTE; 
if  frac( SYMBOLCOUNT  /  4)  =  0.75  then 
TEMPBYTE  :=  (BITS  shl  2)  or  TEMPBYTE; 
if  (frac( SYMBOLCOUNT  /  4)  =  0.0)  then 
begin 

TEMPBYTE  :=  BITS  or  TEMPBYTE; 

TEMPCHAR  :=  chr( TEMPBYTE); 
write( OUTFILE, TEMPCHAR); 

TEMPBYTE:  =0; 
end;  (*if  frac*) 
end;  (*for  I*) 

TEMPCHAR  :=  chr( TEMPBYTE); 

write( OUTFILE, TEMPCHAR);  (*puts  two  0  at  end  of  each  baud*) 
close( OUTFILE); 

XREALa[  k2]  :  =0; 

XIMAGa[ k2] : =0; 
end;  (*DiffDecode*) 


procedure  GraphData(XR,  XI:  TNVectorPtr); 

(*GraphData  graphs  complex  arrays  XR  and  XI.  Two  scales  are  used: Full 
scale,  0  to  kx/2  tones  and  Zoom  scale,  kl  to  k2*) 
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(*loop  indices*) 


var 

grDriver,grMode, ErrCode, 

I,J, 
kl,k2, 

MXDIV2, 

COLOR, 

X, XZ, 

Y, YZ  : integer; 

begin 

kl  :=  round(kx  *  67.0  /  256.0); 
k2  :=  round(kx  *  83.0  /  256.0); 
grDriver  :  =  EGA; 
grMode  : =  EGAHi; 

initgraph(grDriver,grMode, 1 c:  TP4  GRAPHICS'); 

ErrCode  :=  GraphResult; 
if  ErrCode  <>  grOK  then 
begin 

writeln( 'Graphics  error:  ',  GraphErrorMsg(ErrCode)); 
writeln( 'Program  aborted...'); 

Halt(l); 

end;  (*if  ErrCode*) 

SetBKCo lor (black); 

SetColor(white); 

(*axis  cont ruction*) 

line(60,300,600,300);  (*x-axis  for  full  spectrum*) 

line( 68, 140,600, 140);  (*x-axis  for  zoom  spectrum*) 

line(332,300,332,308);  (*x-axis  scale  marks*) 

line(392,300,392,308);  (*x-axis  scale  marks*) 

line(60,300,60,308); 

line(572,300,572,308); 

line(96, 140,332,300); 

line( 5 76, 140,392,300); 

MXDIV2  :=  round(kx  /  2); 

for  I  :=  0  to  NXDIV2  do 
begin 

if  (XRa[ I]  >=  0)  and  (XIa[I]  >  0)  then 
SetColor( lightmagenta); 
if  (XRa[ I]  <  0)  and  (XIa[I]  >  0)  then 
SetColor( lightgreen); 
if  (XRa[I]  <  0)  and  (XIa[I]  <=  0)  then 
SetColor( lightred); 

if  (XRa[I]  >=  0)  and  (XIa[I]  <=  0)  then 
SetColor( yellow); 

(*zoom  spectrum*) 

if  (I  >  kl)  and  (I  <=  k2)  then 
begin 

XZ  :=  round( 2*( ( I  *  4096.0  /  kx)  -  1040)); 

YZ  :=  140  -  round(sqrt(sqr(XRA[ I] )  +  sqr(XlA[ I] ))); 
line(XZ, 140,XZ,YZ); 
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end; 


(-♦full  spectrum*) 

X  :=  round( 2*( I/MXDIV2  *  256  +  30)); 

Y  :=  300  -  round(sqrt(sqr(XRA[ I]  )  +  sqr(XlA[  I]  ))); 
line(X,300,X,Y); 

end; 

(*graph  labels*) 

setcolor( lightgray); 
outtextXY(60,320,'0'); 
outtextXY(331,320, *  kl * ) ; 
outtextXY( 39 1,320, ' k2 1 ); 
outtextXY(565,320, 'kx/21 ); 
outtextXY( 150,330. 1  frequency — >' ); 
outtextXY(20,250, 'amplitude' ); 
readln; 

outtextXY(300, 330, ' Press  return  to  continue.'); 

readln; 

CloseGraph; 
end;  (*GraphData*) 

(* . *) 

(* . *) 

begin  (*main  body*) 

INVERSE  :=  false;  (*sets  forward  FFT*) 

ERROR  :=  0; 

NUMPTS  :=  kx; 

kl  :=  round(kx  *  67.0  /  256.0  +  1); 
k2  :=  round(kx  *  83.0  /  256.0); 

repeat 

writeln( ' Enter  1  to  sample  data'); 
writeln( '  2  to  read  data  from  an  ASCII  file'); 

read ln( ANSWER); 
until  ANSWER  in  [ ' 1* , ' 2'] ; 

if  ANSWER  =  1 1 '  then 

(*AcquireData  samples  input  analog  signal*) 

AcquireData 

else 

GetData; 

(*GetData  reads  a  file  of  discrete  sample  produced  in  the  transmitter. 
GetData  is  used  for  testing  Graphics  and  Decoding  procedures  without 
the  need  for  external  hardware  and  Data  aquisition  board 
initialization*) 

writeln( 'Computing  FFT1); 

Comp lexFFT  ( NUMPTS , INVERSE , XRE AL , XIMAG , ERROR ) ; 

DiffDecode; 

(*  writeln( 'Error  =  ERROR,'  hit  the  enter  key' ); readln; *) 
GraphData(XREAL,  XIMAG); 
end. 
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APPENDIX  J.  RECMES 


program  RECEIVER; 

(*Acquires  the  signal.  Stores  it  in  a  memory  buffer.  Differential  decodes 
between  tones.  Maximum  number  of  bauds  are  received.  The  number  of 
bauds  processed  is  a  user  input*) 

uses  Graph,  Crt,  tp4dl6; 


(*$I-*) 

(*$R-*) 


const  MAX_BUFFER  =  65500; 
type 

TNvector  =  array-0. . 4096-  of  real;  (*TYPE  for  real  and  imaginary  data 

for  FFT  routing*) 

TNvectorPtr  =  ATNvector;  (*Po inter  for  FFT  data  array  which 

allows  dynamic  allocation  of  memory*) 

var 

INVERSE  : boolean; 

X.  AL,  XIMAG  : TNvectorPtr; 

ERROR,  TEMPBYTE: byte; 

J, 

kl,k2,kx, ANSWER, 

ERR_C0DE ,  BAUDCOUNT, 

SYMBOLCOUNT, 


NUMBAUDS, 

MAXNUMBAUDS  :  integer; 

MAGNITUDE, 

PHASE  : real; 

DATAVECTOR  :  a  integer; 

DMAPOINTER  : pointer; 

TESTFILE, 

OUTFILE  :  TEXT; 

(*$I  FFT87B2. INC*) 

(*$I  COMPFFT.  INC*) 

(* . - . *) 

_ _ _ _ _ _ _ *) 


procedure  PacketSetUp; 


begin 

repeat 

clrscr; 

if  kx  <  0  then  writeln('TRY  AGAIN1); 
writeln( 'Enter  baud  size  ' ); 
read ln( ANSWER); 
case  ANSWER  of 
256;  kx:  =  256; 

512:  kx:  =  512; 

1024: kx: =1024; 


66 


2048: kx: =2048; 

4096: kx: =4096; 
end;  (*case*) 
if  kx  =  0  then  kx  :=  -1; 
until  kx  >  0; 

MAXNUMBAUDS  :=  trunc((MAX_BUFFER/2)/kx); 

repeat 

writeln; 

writeln( 'Enter  number  of  ',kx,'  bauds  to  process. 

•  MAXNUMBAUDS,'  is  the  maximum.'); 
r eadlnC  NUMBAUDS ) ; 

until  NUMBAUDS  in  [ 1. . MAXNUMBAUDS] ; 

kl  :=  round(kx  *  67.0  /  256.0  +  1); 
k2  :=  round(kx  *  83.0  /  256.0); 
end;  (*PacketSetUp*) 

(* . 

(* . 

procedure  AcquireData; 

(*AcquireData  initializes  Metrobyte  DASH-16F  data  acquisition  board, 
using  TT00LS  procedure  D16_int  and  D16_ainm.  Data  transfer  is 
controlled  by  the  DMA  controller  and  initialized  by  Dl6_ainm  and 
disabled  by  Dl6_dma_int_disable.  TT00LS  procedures  are  external 
procedures  included  by  'uses'  tp4dl6.*) 

var 

RATE:  real; 

I , CNT_NUM ,  MODE,  CYCLE,  TRIGGER, 

BASE_ADR ,  INT_LEVEL,  DMA_LEVEL, 

*  B0ARD_NUM ,  CHANLO, 

0P_TYPE ,  STATUS,  NEXT.CNT,  ERR_C0DE_S  :  integer; 

►  begin 

B0ARD_NUM  :=  0;  INT_LEVEL  :=  7;  DMA.LEVEL  :=  1; 

BASE_ADR  :=  $300; 

D16_init( B0ARD_NUM , BASE_ADR , INT_LEVEL , DMA_LEVEL , ERR_C0DE ) ; 

CHANLO  :=  0; 

CYCLE:  =0;  (*0-one  sweep  of  the  DMA  1-autoinitialize*) 

TRIGGER: =0;  (*0  -  external  1  -  internal*) 

CNT_NUM:=trunc(MAX_BUFFER  /  2);  (*//  of  samples*) 

RATE  :=  10000. 0; (*used  for  internal  trigger*) 

MODE  :=  2;  (*DMA  mode*) 
writeln( 'Ready  to  acquire'); 

D16_ainm( B0ARD_NUM , CHANLO , MODE , CYCLE , TRIGGER , CNT_NUM , 

RATE ,  DATAVECTORa , ERR_C0DE ) ; 

STATUS  :=  11; 

(*status  indicates  the  progress  of  acquisition.  When  all 
samples  have  been  acquxred  status=0*') 
repeat 

D16_dma_int_status( B0ARD_NUM , 0P_TYPE , STATUS ,NEXT_CNT,  ERR_C0DE_S); 


•*) 

■*) 
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until  STATUS  =  0; 
writeln( 'Data  received'); 

if  ERR_C0DE  <>  0  then 

D 1 6_pr int_er  ror ( ERR_CODE ) ; 

D 1 6_dma_int_dis  ab 1 e ( B0ARD_NUM , ERR_C0DE ) ; 


end; (*Acquire*) 

(* . *) 

(* . . *) 

procedure  ConvertData; 

(*ConvertData  seperates  channel  and  acquired  data.  CHAN_DATA  is  not  used. 
Acquired  data  is  stored  in  XREAL.  *) 

var 

AD_DATA:  array[ 0. . 4096]  of  integer; 

I, CHAN  DATA, ERR  CODE, 

SEGMENTPART , OFFSETPART:  integer; 

NEWDATAVECTOR  : a integer; 

TEMPPOINTER  : pointer; 

begin 

writeln( ' Processing  baud  ' ,BAUDC0UNT); 
fillchar(XREALA ,si2eof  (XREALa) ,0); 
fillchar(XIHAGA ,sizeof  (XIMAGa) ,0); 

SEGMENTPART:  =s  eg( DATAVECTORa ) ; 

OFFSETPART: =ofs(DATAVECT0RA)  +  2  *  kx  *  (BAUDCOUNT  -  1); 

TEMPPOINTER:  =ptr( SEGMENTPART, OFFSETPART); 

NEWDATAVECTOR  :=  TEMPPOINTER; 

dl6_convert_data( 2047 , kx , NEWDATAVECTORa , AD_DATA[ 0] , 

CHAN_DATA , 0 , ERR_C0DE ) ; 

for  I:=  0  to  (kx  -  1)  do 
begin 

XREALa [ I]  : =  AD_DATA[  I] /5; 
end; 

end; (*ConvertData*) 

(* . *) 

procedure  DiffDecode; 

(*DiffDecode  differentially  decodes  complex  frequency  domain  arrays  XREAL 
and  XIMAG.  Four  decoded  symbols  are  recombined  into  a  byte  and 
transferred  to  file  BYTESOUT.DAT.*) 


var 

I  : integer; 

TEMPREAL,TEMPIMAG  :  real; 

BITS  : byte; 

TEMPCHAR  : char; 


begin 

for  I  :=  kl  to  (k2  -  1)  do 
begin 

(*Complex  multiply  two  adjacent  tones,  I  and  the  complex  conjugate  of  1+1. 
This  will  give  the  phase  difference  between  the  two  tones.  The  answer 
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is  in  rectangular  notation*) 

TEMPREAL:  =XREALa[ I]  *  XREALa[  1+1]  +  XIMAGa[  I]  *  XIMAGa[ 1+1] ; 
TEMPIMAG:  =XREALa[ I]  *  XIMAGa[ 1+1]  -  XREALa[  1+1]  *XIMAGa[IJ; 

(*Complex  multiply  (TEMPREAL  +  j  TEMPIMAGE)  and  (Tfj).  This  rotates  the 
differential  vector  pi/4  radians.  XREAL  and  XIMAG  are  used  to  store 
the  results.  This  eliminate  the  original  data*) 

XREALa[ I] :=( TEMPREAL  -  TEMPIMAG)/80; 

XIMAGa[ I] : =(TEMPREAL  +  TEMPIMAG)/80; 


(*decode*) 

if  (XREALa[I]  >=  0)  and  (XIMAGa[I]  >  0)  then 
BITS: =$00; 

if  (XREALa[I]  <  0)  and  (XIMAGa[I]  >  0)  then 
BITS:  =$01; 

if  (XREALa[I]  <  0)  and  (XIMAGa[I]  <=  0)  then 
BITS: =$03; 

if  (XREALa[I]  >=  0)  and  (XIMAGa[ I]  <=  0)  then 
BITS: =$02; 

SYMB0LC0UNT  :  =  SYMB0LC0UNT  +  1; 

( * f i 1 1  TEMPBYTE  with  four  symbols*) 
if  f rac( SYMBOLCOUNT  /  4)  =  0.25  then 
TEMPBYTE  :=  (BITS  shl  6); 
if  f rac( SYMB0LC0UNT  /  4)  =  0. 5  then 
TEMPBYTE  :=  (BITS  shl  4)  or  TEMPBYTE; 
if  f rac( SYMBOLCOUNT  /  4)  =  0.75  then 
TEMPBYTE  :=  (BITS  shl  2)  or  TEMPBYTE; 
if  (frac( SYMBOLCOUNT  /  4)  =  0.0)  then 
begin 

TEMPBYTE  :=  BITS  or  TEMPBYTE; 

TEMPCHAR  :=  chr( TEMPBYTE); 
wr it  e( OUTF I LE , TEMPCHAR ) ; 

TEMPBYTE: =0; 
end;  (*if  frac*) 
end;  (*for  I*) 


XREALa[ k2] : =1. 0; 

XIMAGa[  k2] :  =1.  0; 

end;  (*DiffDecode*) 

(* . *) 

procedure  Showmessage; 

(*Showmessage  reads  in  decoded  message*) 


var 

NEXTCHAR: char; 

begin 

writeln; 

writeln('The  message  transmitted  is..  '); 
assign(OUTFILE, 'MESSAGE.  DAT' ); 
reset(OUTFILE); 
while  not  EOF(OUTFILE)  do 
begin 

while  not  EOLN(OUTFILE)  do 
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begin 

r  e  ad ( OUTF I LE , NEXTCHAR ) ; 
write(NEXTCHAR) ; 
end; (*while  not  EOLN*) 
readln(OUTFILE); 
writeln; 

end; (*while  not  EOF*) 
close(OUTFILE); 
end; (*Showmessage*) 

(* . *) 

begin  (*main  body*) 

GetDMABuf f er (MAX_BUFFER , DMAPO INTER , ERR_C0DE); 

DATAVECTOR  :=  DMAPOINTER;  (*This  statement  assigns  a  generic  pointer 

to  a  variable  of  a  specific  pointer  type, 
i. e.  Ainteger,  so  that  the  pointer  can  be 
passed  to  the  dl6_ainm  routine.*) 


new(XREAL); 

new(XIMAG); 


INVERSE  :=  false;  (*sets  forward  FFT*) 

ERROR  :=  0; 
kx:  =0; 

PacketSetUp; 

(*received  message  output  file*) 
assign(OUTFILE, 'MESSAGE.  DAT' ); 
rewrite(OUTFILE); 

(*RECSTAT  file  shows  the  computed  real  and  imaginary  values  of  the 
received  signal  and  its  decoded  representation*) 
assign(TESTFILE, 'RECSTAT.  DAT' ); 
rewrite(TESTFILE); 

(*  writeln(TESTFILE, ' Baudlength  is',kx:5);*) 

SYMB0LC0UNT:  =0; 

TEMPBYTE: =0; 

AcquireData;  (*AcquireData  samples  input  analog  signal*) 

for  BAUDCOUNT  :  =  1  to  NUMBAUDS  do 
begin 

ConvertData; 

Compl exFFT  ( kx , INVERSE , XREAL , XIMAG , ERROR ) ; 

(*  writeln(TESTFILE , ' Baud ' :  11 Tone ' :  5 . ' Real ' :  10, 

' Imaginary' :  10, 'Mag' : 10, 'Phase' : 8); *) 
for  J: =kl  to  k2  do 
begin 

MAGNITUDE:  =20*ln(sqrt( sqr(XREALA[ J]  )+ 
sqr( XIMAGa [  J] ) ) )  /  ln(10); 

if  (XREALa[J]  >=  0)  then 
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PHASE: =arctan(XIMAGA[ J]  /  XREALa[J]); 
if  (XREALa[J]  <  0)  and  (XIMAGa[J]  >  0)  then 
PHASE: =Pi  +  arctan(XIMAGA[ J]  /  XREALa[J]); 
if  (XREALa[J]  <  0)  and  (XIMAGa[J]  <=  0)  then 
PHASE:  =arctan(XIMAGA[  J]  /  XREALa[J])  -  Pi; 

writeln(TESTFILE, 'R' ,BAUDC0UNT:  4, J:  5 , 

XREALa[  JJ:  10:4,XIMAGa[  J] :  10:4, 
MAGNITUDE:  10:  4, PHASE:  8:  4); 

end; 

DiffDecode; 

for  J: =kl  to  k2  do 
begin 

MAGNITUDE:  =sqrt(sqr(XREALA[J]  )+sqr(XIMAGA[  J]  )); 
write ln(TESTFILEs 'D' , BAUDCOUNT:  4 , J:  5 , 

XREALa[ J) /MAGNITUDE:  10:  4, 
XIMAGa[ J] /MAGNITUDE:  10:  4); 

end; 

end; 

writeln(OUTFILE , chr( TEMP BYTE)) ; 

close(OUTFILE); 

close(TESTFILE); 

dispose(XREAL); 

dispose(XIMAG); 

Fr eeDMABuf  f er ( MAX_BUFFER , DMAPOINTER , ERR_C0DE ) ; 

Showmessage; 

(*  writeln(YError  =  'jERROR,'  hit  the  enter  key' );  readln; *) 
end. 
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APPENDIX  K.  RESPONSE 


program  RESPONSE; 

(*RESP0NSE  outputs  the  channels  phase  response,  by  subtracting 
the  transmitted  phase  from  the  received  phase  for  each  tone. 
Prior  to  running,  set  NUMTONES  and  K1MINUS1  for  the  proper 
baud  parameters  then  compile  and  run  the  program*) 

var 

MAG, PHASE, 

XREAL,XIMAG  : array[ 1. . 256]  of  real; 

K1MINUS1, 

NUMTONES, K, 

I,BAUD,Kx,Kr  : integer; 

XMITPHASE, 

XR,XI  : real; 

CODE  : char; 

XMITDATA, 

RECDATA , 

CHANNELDATA  : text; 
begin 

as sign( XMITDATA, 'xmitdatl. 024' ); 
reset (XMITDATA); 

ass ign( RECDATA, ' recstat. dat' ); 
reset(RECDATA); 

ass ign( CHANNELDATA, 'respl024. dat' ); 
rewrite( CHANNELDATA) ; 

(*NUMT0NES  and  K1MINUS1  should  be  set  for  the  same  baud  parameters 
as  xmitdat*. ***  and  recstat.dat*) 

NUMTONES: =64; 

K1MINUS1: =268; 

(*get  transmitted  phase  and  received  phase*) 
while  not  E0F( XMITDATA) do 
begin 

for  I: =1  to  NUMTONES  do 
begin 

readln( XMITDATA, BAUD,Kx,XREAL[ I] ,XIMAG[ i] ); 
readln(RECDATA, CODE, BAUD, Kr,XR, XI ,MAG[ I] ,PHASE[I] ); 
end; 

(*ensure  tones  match*) 
if  Kx  <>  Kr  then 
begin 

writeln('Kx  <>  Kr'); 
readln; 
end; 

(*convert  transmitted  symbol  into  phase*) 
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for  I: =1  to  NUMTONES  do 
begin 

if  (XREAL[I]>=  0.0)  and  (XIMAG[I]>=  0.0)  then 
XMITPHASE;  =  pi/4; 

i£  (XREAL[  I]  <  0.0)  and  (XIMAG[i]>=  0.0)  then 
XMITPHASE:  =  3*pi/4; 

if  (XREAL[  I]  >=  0.0)  and  (XIMAG[I]<  0.0)  then 
XMITPHASE:  =  -pi/4; 

if  (XREAL[  I]  <  0.0)  and  (XIMAG[I]<  0.0)  then 
XMITPHASE:  =  -3*pi/4; 

(*subtract  transmitted  phase  from  received  phase*) 

PHASE[ I] :  =PHASE[ I]  -  XMITPHASE; 

(*phase  differences  may  jump  be  2*pi  radiar  This  helps 
but  does  not  work  in  every  case.  Output  file  may  have 
to  be  edited  to  2*pi  jumps*) 

0  if  PHASE [ I]  >  2*pi  then 

PHASE[ i] :  =PHASE[  I]  -  2*pi; 
if  PHASE[  ij  <  -2  then 

PHASE[ i] :  =PHASE[ i]  +  2*pi: 

(*actual  tone  numbers*) 

K:=I  +  K1MINUS1; 

writeln(CHANNELDATA,BAUD:  4,K:  4,MAG[  I] :  10:  4,PHASE[  .  10:  4); 
readin(RECDATA); 
end; 

end; 

clo-,e(XMITDATA); 
clc^eCRECDATA); 
c los  e( CHANNELDATA) ; 

erd. 
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APPENDIX  L.  STATISTICS 


program  Statistics; 

(*to  run  this  program  the  transmitter  program  Xmitmes  must  be  run  using 
MESSAGE.DAT.  The  file  produced,  XMITDAT.DAT,  must  be  transferred  to  the 
receiver  computer.  Then  transfer  the  data  between  computers  using 
Xmitmes  and  Receive.  Receive  produces  the  file  RECSTAT.  DAT.  This 
program  uses  XMITDAT.DAT  and  RECSTAT.DAT  to  compute  the  mean  and 
variance  for  each  quadrant,  and  the  SNR  for  each  quadrant,  baud  and 
the  overall  SNR  for  the  packet.*) 

uses  CRT; 

type 

STATARRAY  =  array[ 1. . 125]  of  real; 


var 

Kx , K , NUMBAUDS , 

I, J, BAUD, 

TONEX, TONER, 

Q1,Q2,Q3,Q4, 

BITERR0R1 , BITERR0R2 , BITERR0R3 , BITERR0R4  : integer; 

TEMPR, TEMPI, 

SNRINDB,SNRIN, 

SNROUTDB , SNROUT, 

MEANRjMEANI ,VARR,VARI: real; 

CODE  : char; 

XD,RD,SD  : text; 

XR,XI ,RR,RI  :array[l.  .  16,1.  .  128]  of  real; 

BTM , BTV , BTQ  : array[ 1. . 16]  of  real; 

SR1,SR2,SR3,SR4, 

SI1,SI2,SI3,SI4  :  STATARRAY; 

STATMATRIC  :  array[  1. .  16,1. .  4,1. .  5]  of  real; 


(* . *) 

procedure  Stat(NPTS: integer; 

X:  STATARRAY; 

var  XMEAN ,  XVAR  :real); 

(*input  NPTS:  number  of  points  to  be  used  to  compute  the  mean 
X: the  array  of  data 

output  XMEAN:  the  mean  of  the  data  in  array  X 

XVAR: the  variance  of  the  data  in  array  X*) 


SUM: real; 

N  : integer; 
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begin 

(*Compute  mean*) 

SUM:  =  0.0; 

for  N:=  1  to  NPTS  do 
SUM:  =  SUM  +  X[N] ; 

if  NPTS  >  0  then 

XMEAN: =  SUM  /  NPTS 
else 

XMEAN:  =1.  Oe- 10; 

( ’’'Compute  variance*) 

SUM:  =  0.0; 

for  N: =1  to  NPTS  do 

SUM:  =  SUM  +  sqr(X[ N]  -  XMEAN); 

if  NPTS  >  1  then 

XVAR:  =  SUM  /  (NPTS  -  1) 

elSXVAR:=1.0e-10; 

end;  (*Stat*) 

(••'•' . 

begin  (*ma.in  body*) 
repeat 

write('Enter  baud  size.  '); 
readln(Kx); 

until  (Kx  mod  256  =  0)  and  (Kx  >=  256)  and  (Kx  <=  4096); 

writeC  Enter  input  SNR  in  DB  ' ); 
readln(SNRINDB); 

SNRIN  :=  exp(SNRINDB  /  10  *  ln(10)); 

K  :=  trunc(Kx  /  16); 

assign(XD, 'XMITDAT.  DAT' ); 
reset(XD); 

assign(RD, ' RECSTAT.  DAT1 ); 
reset(RD); 
case  Kx  of 

256: begin  assign(SD, ' STAT256.  DAT'  );  NUMBAUDS 
512: begin  assign(SD, ' STAT512.  DAT' );  NUMBAUDS 
1024: begin  assign(SD, ' STAT1024.  DAT’ );  NUMBAUDS 
2048: begin  assign(SD, ' STAT2043.  DAT' );  NUMBAUDS 
4096: begin  assign(SD, ' STAT4096.  DAT' );  NUMBAUDS 
end; (*case  Kx*) 
rewrite(SD); 

(*read  in  transmitted  and  received  data*) 
for  I: =  1  to  NUMBAUDS  do 
begin 

for  J:  =  1  to  K  do 
readln(RD); 
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for  J: =  1  to  K  do 
begin 

readln(XD,  BAUD,  TONEX,  XR[I,J],  XI[  I,J]  ); 
readln(RD,  CODE,  BAUD,  TONER,  RR[I,J],  RI[I,J])j 
if  TONEX  <>  TONER  then 
begin 
clrscr; 

writeln( 'XMITDAT.  DAT  and  RECSTAT.DAT  are  in  error.'); 
Halt; 
end; 

end; (*for  J*) 
end; (*for  I*) 
close(XD); 
close(RD); 

(*decode  transmitted  data*) 
for  I:-  1  to  NUMBAUDS  do 
for  J:=  1  to  (K  -  1)  do 
begin 

TEMPR:  =  XR[  I ,  J]  *XR[I,J+1]  +XI[I,J]  *XI[I,J+1]; 

TEMPI: =  XR[  I , J]  *  XI[I,J+1]  -  XR[I,J+1]  *  XI[  I,J] ; 

XR[ I, J]:=  TEMPR  -  TEMPI; 

XI[ I , J] :  =  TEMPR  +  TEMPI; 
end; (*for  J*) 

(*store  statistics  for  each  quadrant  and  baud  compute  statistics*) 
BITERR0R1: =0;  BITERR0R2: =0;  BITERR0R3: =0;  BITERR0R4:  =0; 
for  I:=  1  to  NUMBAUDS  do 
begin 

Ql:  =0;  Q2:  =0;  Q3:=0;  Q4:  =0; 
for  J:=  1  to  (K  -  1)  do 
begin 

if  (XR[  I ,  J]  >=  0.0)  and  (XI[  I, J]  >  0.0)  then 
begin 

Q1:=Q1  +  1; 

SR1[ Ql] :  =  RR[  I,J] ; 

SI1[Q1]:=  RI[  I,J] ; 
if  RR[ I, J]  <  0. 0  then 

BITERR0R1:  =  BITERROR 1  +  1; 
if  RI[  I ,  J]  <=  0.0  then 

BITERR0R1:  =  BITERR0R1  +  1; 

end; 

if  (XR(  I , J]  <  0.0)  and  (XI[I,JJ  >  0.0)  then 
begin 

Q2:=Q2  +  1; 

SR2[ Q2] :  =  RR[I,J]; 

SI2[Q2]  :=  RI[  I, J] ; 
if  RR[  I, J]  >=  0.  0  then 

BITERR0R2:  =  BITERR0R2  +  1; 
if  RI[I,J]  <=  0.0  then 

BITERR0R2:  =  BITERR0R2  +  1; 

end; 

if  (XR[  I ,  J]  <  0.0)  and  (XI[I  ,Jj  <=  0.0)  then 
beein 

~Q3:  =Q3  +  1; 

SR3[Q3]:=  RR[  I,J] ; 
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SI3[Q3]:=  RI[  I, J] ; 
if  RR[I,J]  >=  0.0  then 
B I TERROR 3: =  BITERR0R3  + 
if  RI[ I ,  J]  >0.0  then 

BITERR0R3: =  BITERROR3  + 

end; 

if  (XR[  I, J]  >=  0.0)  and  (XI[  I , J] 
begin 

Q4:=Q4  +  1; 

SR4[Q4]:=  RR[  I ,  J]  ; 

SI4[  Q4] :  =  RI[  I ,  J] ; 
if  RR[  I, J]  <0.0  then 

BITERR0R4:  =  BITERR0R4  + 
if  RI[ I, J]  >  0.  0  then 

BITERR0R4:  =  BITERR0R4  + 

end; 

end;  (*for  J*) 


l; 

i; 

<=  0. 0)  then 


i; 

i; 


(*STATMATRIC  is  a  3  dimensional  array.  The  first  index  is  the  baud 
number.  The  second  is  the  quadrant.  The  third  is  the  statistics 
for  each  quadrant. 

.°T^TMATRIC[  I, J,  1]  =  the  mean  of  the  real  parts 

AT^Tcf  I, J,2]  =  the  variance  of  the  real  parts 
~V  ,VuC[  I,J,3]  =  the  mean  of  the  imaginary  parts 
°V.  MRIC[I,J,4}  =  the  variance  of  the  imaginary  parts 

’  .jATRIC[ I, J,5]  =  the  number  of  points  transmitted  in  a  given 
quadrant*) 


BTM[  I]  :  =0.  0; 

BTV[  I]  :  =0.  0; 

BTQ[  I]  :  =0.  0; 
if  Q1  >  1  then 
begin 

S t  a t ( Q 1 , SR 1 , MEANR , VARR ) ; 

Stat( Q1 ,  SI  1 , MEANI , VARI ) ; 

STATMATRIC[ 1,1,1]:  =MEANR; 

STAiMATRIC[ 1,1,2]  :=VARR; 

STATMATRiq 1,1,3]: =MEANI; 

STATMATRIC[ 1,1,4]:  =VARI; 

BTM[ I] :  =BTH[ I]  +  Ql*( abs ( MEANR) +abs ( MEANI ) ) ; 
BTV[ I]  : =BTV[ I]  +  (Q1-1)*(VARR+VARI); 

BTQ[  I] :  =BTQ[  I]  +  2*Q1; 
end; 

STATMATRIC[ 1,1,5]  :=Q1; 

if  Q2  >  1  then 
begin 

Stat(Q2,SR2, MEANR , VARR ) ; 

Stat(Q2,SI2, MEANI, VARI); 

STATMATRIC( 1,2,1]:  =MEANR; 

STATMATRIC[ 1,2,2]: =VARR; 

STATMATRICf 1,2,3]:  =MEANI; 

STATMATRIC[ 1,2,4]:  =VARI; 

BT:1[  I] :  =BTM(  I]  -f  Q2*(abs(MEANR)+abs(MEANI)); 
BTV[ I] :  =BTV[  I]  +  (Q2-1)*(VARR+VARI); 

BTQ[  I]  :  =BTQ[  I]  +  2*Q2; 
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end; 

STATMATRIC[ 1,2,5): =Q2; 

if  Q3  >  1  then 
begin 

S  t  at ( Q3 , SR  3 , ME ANR , VARR ) ; 

Stat( Q3 , SI3 , MEANI , VARI ) ; 

STATMATRiq 1,3,1]:  =MEANR; 

STATMATRiq 1,3,2]:  =VARR; 

STATMATRiq 1,3,3]:  =MEANI; 

STATMATRIC[ 1,3,4]:  =VARI; 

BTM[ I] : =BTM[ I]  +  Q3*( abs ( MEANR)+abs (MEANI ) ) ; 
BTV[ I] : =BTV[ I]  +  (Q3-1)*(VARR+VARI); 

BTQ[ I] : =BTQ[ I]  +  2*Q3; 
end; 

STATMATRiq  1,3,5]  :=Q3; 

if  Q4  >  1  then 
begin 

Stat( Q4 , SR4 , MEANR , VARR) ; 

Stat(Q4,SI4, MEANI, VARI); 

STATMATRiq  1,4,1]:  =H£ANR; 

STATMATRiq 1,4,2]:  =VARR; 

STATMATRIC[ 1,4,3]:  =MEANI; 

STATMATRIC[ 1,4,4] : =VARI; 

BTM[ I]  :=BTM[  I]  +  Q4*( abs ( MEANR) +abs ( MEANI ) ) ; 
BTV[ I] : =BTV[ I]  +  (Q4-1)*(VARR+VARI); 

BTQ[ I]  : =BTQ[ I]  +  2*Q4; 
end; 

STATMATRiq  1,4,5]  :  =Q4; 
end; (*for  I*) 


(*output  statistics  to  a  file*) 

writeln(SD, '***  SNRIN=  ’ ,SNRIN: 10: 7, '  or  ’ ,SNRINDB:  10: 7, '  DB  ***' ); 
writeln(SD); 

writeln(SD,'  Baud  length  is  1 ,Kx); 

SR2[  1] :  =0.  0; 

SR2[  2] :  =0.  0; 

SR2[  3] :  =0.  0; 

for  I:=  1  to  NUMBAUDS  do 
begin 

writeln(SD);  writeln(SD); writeln(SD); writeln(SD); 
writeln(SD, 'DFT  statistics  for  baud  #:  1 ,1: 4); writeln(SD); 
for  J:=  1  to  4  do 
begin 

(*real  output  data*) 

writeln(SD); writeln(SD); 

writeln(SD, ' Given  the  phase  is  in  quadrant  ' ,J); 
writeln(SD); 

if  STATMATRiq I, J, 5]  >  1  then 
begin 

writeln(SD, 'Number  of  points  =  ' :  35 , 

STATMATRiq  I, J, 5] :  10:  5); 
writeln(SD, 'Mean  of  real  parts  =  ' : 35 , 

STATMATRiq  I ,  J,  1]  :  10:  5); 
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writeln(SD, 'Variance  of  real  parts  =  ’ : 35 , 
STATMATRIC[I,J,2] :  10:5); 

SNROUT:  =s qr ( STATMATR I C [ I , J, 1] )  /  STATMATRIC[ I , J,2J ; 
SNROUTDB:  =  10.0  *  ln( SNROUT)  /  ln(10. 0); 
writeln(SD, ' SNROUT  of  real  part  =  ’:35, 

SNROUT:  10:5,'  =  ',  SNROUTDB: 10: 5, '  DB'); 
writeln(SD, ’ SNR0UT/SNRIN  of  real  part  =  ’:35, 

►  (SNROUT  /  SNRIN): 10:5); 

(*imaginary  output  data*) 
v  writeln(SD); 

writeln(SD, 'Mean  of  imag  parts  =  ' : 35 , 

STATMATRIC[ I ,  J,  3] :  10:5); 
writeln(SD, 'Variance  of  imag  parts  =  ':35, 
STATMATRIC[ I , J,4] :  10:  5); 

SNROUT: =  sqr(STATMATRIC[ I , J,3] )  /  STATMATRIC[ I , J,4] ; 
SNROUTDB:  =  10  *  ln( SNROUT)  /  ln(10.0); 
writeln(SD, ' SNROUT  of  imag  part  =  ':35, 

SNROUT:  10:  5 , '  =  SNROUTDB: 10:  5, '  DB* ); 
writelnCSD,' SNROUT/ SNRIN  of  imag  part  =  ’ :  35 , 

(SNROUT  /  SNRIN): 10:5); 

end  (*if  Statmatric*) 

else 

begin 

writeln(SP, 'Quadrant  has  less  then  2  statistic 
points  '); 

writeln(SD, ' and  is  not  used  in  the  overall  statics.'); 
end;  (*else*) 

»  end; (*for  J*) 

writeln(SD); writeln(SD);  writeln(SD); 

writeln(SD,  'Overall  (Real  +  Imag)  statistics  for  Baud  #:  '  ,  1:4); 

SNROUT:  =  sqr(BTM[ I]  /BTQ[ I] )  /  (BTV[ I] /(BTQ[ I]  -1)); 

SNROUTDB: =  10  *  ln( SNROUT)  /  ln(10.  0); 

writeln(SD, 'Overall  Baud  SNROUT  =  ’ :  35 , SNROUT:  10:  5 , ’  = 

SNROUTDB:  10:5,’  DB1); 

writeln(SD, ’Overall  baud  gain  (SNROUT  /  SNRIN)  =  ’ : 35 , 

(SNROUT/ SNR IN):  10:5); 

(*store  baud  means  and  variances*) 

SR2(  1] :  =  SR2[  1]  +  BTM[  I] ; 

SR2[  2] :  =  SR2[  2]  +  BTV[  I] ; 

SR2[  3] :  -  SR2[  3]  +  BTQ[  I] ; 
end; (*for  I*) 

writeln(SD); writeln(SD);  writeln(SD);  writeln(SD); 

writeln(SD, ’***  TOTAL  OVER  ALL  1 , NUMBAUDS , 1  BAUDS  ABOVE  BAUD  ', 

Kx, '  ***'); 

writeln(SD);  writeln(SD); 

SNROUT:  =sqr(SR2[l]/SR2[  3])  /  ( SR2[  2]  /  ( SR2[  3]  -  1)); 

SNROUTDB:  =10  *  Inf SNROUT)  /  ln(10. 0); 

writeln(SD, ' Overall  SNROUT  =  ':  35, SNROUT:  10:  5, '  =  ', 

SNROUTDB:  10:5,'  DB'); 
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writelnC 'Overall  SNROUT  =  ': 35. SNROUT:  10:  5, '  = 

SNROUTDB: 10:5,'  DB ' ) ; 

writelnC SD, 'Overall  gain  (SNROUT  /  SNRIN)  =  ' :  35 , 

( SNROUT/ SNRIN): 10:5); 

writelnC SD) ; writelnC SD) ; 

writelnC SD, 'Total  possible  statistical  points  ' , (Kx/8-2)*NUMBAUDS: 3: 0, 
Actual  number  of  statistical  points  ' ,SR2[ 3] : 3: 0); 


writelnC SD); 

writelnC SD, 'BIT  ERRORS': 35); 
writelnC SD, 'Quadrant  1  =  ' : 30,BITERR0R1); 
writelnC SD, ' 2  =  ’ :  30,BITERR0R2); 
writeln(SD, ' 3  =  ' :  30,BITERR0R3); 
writelnC SD, '4  =  ’ :  30,BITERR0R4); 

BITERROR 1: =  BITERROR 1  +  BITERR0R2  +  BITERROR3  +  BITERR0R4; 
writelnC SD, 'Total  Bit  errors  =  * : 35 ,  BITERR0R1); 
writelnC 'Total  Bit  errors  =  ' : 30 ,  BITERR0R1); 
close(SD); 
end. (*main  body*) 
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